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I.  INTRODUCTION  AND  SUMMARY 


This  report  describes  the  results  obtained  during  the  second  year  of 
a three  year  study  of  the  ocean  bottom  interaction  of  underwater  sound. 

The  work  reported  here  is,  for  the  most  part,  still  in  progress;  thus, 
this  report  portrays  an  ongoing  investigation,  not  the  results  of  a com- 
pleted study. 

For  this  reason,  among  others,  this  report  is  focused  on  the  scientific 
problem  of  acoustic  bottom  interaction  itself  and  not  on  tlie  intended 
applications.  Beginning  in  the  third  year  of  this  study  these  results, 
as  well  as  the  results  obtained  earlier,  will  be  brought  to  bear  on 
questions  oriented  to  specific  applications. 

The  results  appearing  herein  include  material  which  has  not  yet 
been  reported  elsewhere,  as  well  as  brief  summaries  of  work  which  have 
already  appeared  in  specialized  technical  reports,  papers,  and  oral 
presentations.  This  appearance  of  original  material,  as  well  as  the 
necessity  for  documenting  certain  aspects  of  this  research  program, 
accounts  for  the  rather  high  level  of  detail  and  completeness  of  this 
report . 

A.  Background:  The  Acoustic  Bottom  Interaction  Problem 

Until  comparatively  recently.  Interest  in  the  Interaction  of  underwater 
sound  with  the  ocean  bottom  was  confined  to  shallow  water,  high  frequency 
bottom  reflection  loss  studies  by  acousticians  and  short  range  deep 
acoustic  probing  of  the  ocean  floor  by  seismologists.  As  the  attention 
of  acousticians  has  been  drawn  toward  low  frequency  (below  1 kHz)  long 
range  propagation  problems,  the  necessity  has  grown  to  better  understand 
and  describe  the  details  of  sound  propagation  in  and  near  the  ocean  floor, 
especially  the  deep  ocean  basins  and  their  defining  slopes  and  ridges.  At 
the  same  time,  the  surge  of  interest  and  effort  in  the  area  of  plate 
tectonics  has  begun  to  make  available  increasingly  detailed  information 
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concerning  the  subbottom  structure  of  these  areas. 


The  acoustic  bottom  interaction  problem  can  be  broadly  defined  as 
quantifying  the  importance  of  the  ocean  bottom  substructure  in  determining 
the  characteristics  of  propagated  sound  fields,  regions  in  the 
water  column  of  comparative  insensitivity  to  bottom  interaction, 
the  exploitability  of  bottom  traveling  energy,  etc.  The  range  of 
such  problems,  considering  the  variability  of  ocean  bottom  substructure, 
acoustic  frequencies,  source-receiver  geometries,  etc.,  is  considerable 
and  it  is  clear  that  one  immediate  problem  is  to  define  a series  of 
specific  questions  which  can  be  usefully  addressed. 

Tiie  simplest  subset  of  problems  which  still  has  a significant 
bearing  on  the  overall  problem  and  importance  to  the  applications  is 
obtained  by  assuming  a horizontally  stratified  subbottom  structure  and 
using  the  associated  plane  wave  reflection  coefficient  (bottom  loss)  as 
an  investigative  tool.  Such  bottom  loss  studies  constitute  an  integral 
part  of  this  research  program.  Questions  concerning  the  depth  to  which 
any  subbottom  information  is  required  (hidden  depth  problem),  the  importance 
of  shear  waves  as  a loss  mechanism,  the  accuracy  with  which  various 
parameters  must  be  specified,  etc.,  may  be  investigated  without  the  addi- 
tional complications  of  specifying  source-receiver  geometry,  (water) 
sound  speed  profile,  and  bathymetry. 

There  are,  of  course,  many  important  questions  which  do  involve 
these  variables  in  an  integral  way,  such  as  those  Involving  isolating 
regions  in  the  water  column  of  relative  insensitivity  to  bottom  interaction, 
the  interplay  between  the  refractive  effects  of  the  water  sound  speed 
structure  and  bottom  interaction,  the  source-receiver  range  dei^endence  of 
the  partitioning  of  propagated  energy  between  bottom  interacting  and 
purely  waterborne  paths,  etc.  These  questions,  as  well  as  others,  can  be 
usefully  addressed,  in  an  approximate  fashion,  by  assuming,  again,  horizon- 
tal stratification  and  Ignoring,  for  the  moment,  the  effects  of  range 
variable  bathymetry  and  sound  speed  structure.  A variety  of  problems  of 
this  type,  some  of  which  are  discussed  in  this  report,  are  under 
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investigation  as  part  of  this  research  program. 


A third  level  of  detail  involves  problems  associated  with  range 
changing  bathymetry,  propagation  over  continental  slopes,  sea  mounts, 
ridges,  etc.,  as  well  as  the  associated  problems  of  interface  roughness. 
These  problems  are  among  the  most  difficult  in  this  entire  area  and  must 
themselves  be  addressed  piecemeal.  It  is  useful,  for  example,  to  ignore 
subbottom  structure  and  employ  a phenomenological  bottom  reflection  loss 
within  a ray  trace  model.  This  can  then  be  used  to  Investigate  questions 
concerning  the  importance  of  bottom  reflectivity  in  multipath  conversion 
process  associated  with  propagation  over  continental  slopes  (e.g.,  slope 
enhancement ) . 

Another  class  of  problems  in  this  area  concerns  the  sensitivity  of 
propagation  characteristics  to  bottom  roughness,  either  water-sediment 
interface  roughness  or,  perhaps  more  Importantly,  subbottom  (basement) 
roughness.  One  possibility  for  investigating  such  questions  could  be 
to  employ  the  results  of  rough  surface  scattering  research  to  estimate  the 
scattering  loss.  Such  an  approach,  if  applicable,  would  in  fact  be  an 
extension  of  the  bottom  loss  studies  discussed  earlier.  A second  approach 
would  be  to  include  the  roughness  directly  in  the  propagation  problem, 
such  as  by  the  addition  of  an  appropriate  bottom  impedance  modification. 
Both  of  these  approaches,  as  well  as  others,  have  been  considered  and 
are  being  piirsued. 

The  preceding  breakdown  of  the  bottom  interaction  problem  has 
involved  a consideration  of  single  sources  and  point  receivers.  Some  of 
the  methods  used  in  these  studies  can  be  extended  to  incliide  consideration 
of  multiple  sources  and  receivers.  The  objectives  of  such  an  extension 
are  twofold,  to  Investigate  the  bottom  interaction  effects  on  more  complex 
acoustic  fields,  such  as  ambient  noise  (or  a combination  of  signal  and 
noise  fields),  as  well  as  to  consider  questions  of  spatial  coherence 
and  how  it  is  affected  by  bottom  substructure,  bathymetry  changes,  etc. 
Both  of  these  topics  constitute  Important  areas  of  concern  and,  although 
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not  considered  in  this  report,  constitute  active  aspects  of  ongoing  and 
future  research. 

B.  Structure  of  the  ARL:UT  Bottom  Interaction  Research  Program 

The  breakdown  of  the  general  bottom  interaction  problem  into  more 
specific  components,  each  designed  for  the  investigation  of  certain 
physical  mechanisms,  was  in  fact  the  most  important  outgrowth  of  the  first 
year  of  this  research  program.  Much  of  the  first  year's  effort  (already 
reported)  went  into  the  planning  and  outlining  of  a series  of  problems 
which  could  be  addressed  in  a practical  way  and  the  sum  total  of  which 
would  constitute  the  solution  of  an  important  and  applicable  segment  of 
the  bottom  interaction  problem. 

During  the  first  year's  effort  the  major  lines  of  investigation 
were  begun  and  preliminary  results  were  obtained,  primarily  via  model 
studies  using  existing  models  obtained  from  other  laboratories.  During 
the  course  of  those  investigations  it  became  clear  that  more  specialized 
calculational  techniques,  appropriate  to  the  then  defined  range  of  problems, 
would  be  required.  Development  of  these  methods,  chiefly  ray  trace  methods, 
a normal  mode  model,  and  a bottom  reflection  loss  model,  was  begun. 

During  the  subsequent  year,  the  period  reported  here,  this  developmental 
effort  was  completed  and  extensive  application  of  these  tools  was  under- 
taken. Continued  use  was  made  of  existing  models,  although  in  some  cases 
considerable  modification  was  required. 

During  the  first  two  years  of  this  program,  the  second  of  which  is 
reported  here,  the  emphasis  has  been  primarily  on  the  scientific  bottom 
interaction  problem  itself,  with  only  general  guidance  providing  the  link 
to  applications.  During  the  third  year  a major  component  of  the  study 
will  be  to  assemble  the  individual  components  to  form  an  integrated 
series  of  recommendations  and  results  oriented  toward  specific  applications. 
With  this  plan  in  mind,  little  attention  is  paid  in  this  report  to  these 
matters. 
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Accompanying  the  theoretical  parameter  studies  using  models  has  been 
an  effort  to  employ  the  results  of  data  analysis  and  interpretation  to 
these  problems.  This  effort  began  dm-ing  the  first  year  of  this  program 
and  has  continued  throughout.  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin  (ARL:OT)  acoustic  data  archives  were  initially 
composed  only  of  data  acquired  from  LRAPP  sponsored  exercises,  and  more 
recently  expanded  to  include  data  from  exercises  sponsored  by  Naval  Electronic 
Systems  Command  (NAVEIEX),  Code  FME-22h,  as  well  as  newly  acquired  LRAPP 
data.  These  data,  together  with  a modeling  capability  specifically  suited 
to  the  investigation  of  bottom  interaction  problems,  constitute  the  primary 
investigative  resoiirces  of  this  program. 


The  processing,  analysis,  and  intrepretation  of  both  existing  data 
and  newly  acquired  data  form  an  integral  element  of  various  ARL:UT  ongoing 
research  programs,  the  spinoff  from  which  has  served  as  an  important 
source  of  guidance  and  direction  for  the  bottom  interaction  program,  sponsored 
by  NAVEIEX,  Code  520.  These  additional  programs  include  ambient  noise 
(NAVELEX,  Code  320,  and  NORDA/LRAPP,  Code  600),  sensor  performance,  RDSS 
(NAVELEX,  Code  FME-124),  multipath  analysis  and  bottom  limited  propagation 
(NAVELEX,  Code  FME-124  and  NORDA/LRAPP,  Code  6oo),  and  in  situ  sediments 
properties  (NORLA,  Code  48o). 

C.  Contents  of  the  Report 

The  technical  portion  of  this  report  is  divided  into  fiVe  additional 
chapters,  the  first  three  of  which.  Chapters  II  through  IV,  are  concerned 
with  the  major  research  areas  of  this  program. 


Chapter  II  considers  a variety  of  bottom  interaction  questions 
addressable  by  using  the  bottom  reflection  loss  as  a measure  of  the  impor- 
tance of  variations  in  subbottom  description.  The  first  section  of  this 
chapter  considers  the  importance  of  the  presence  of  a (basalt)  substrate, 
or  basement,  to  the  bottom  reflection  loss.  In  particular,  the  existence 
of  a solid  substrate  under  the  fluid  sediment  layers  Introduces  the 
possibility  of  a boundary,  or  interface  wave,  called  a Stoneley  wave  which 
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can  act  as  an  additional  loss  mechanism.  The  second  and  third  sections 
of  Chapter  II  consider  the  importance  of  depth  variations  (gradients)  of 
attenuation  and  density  within  sediments.  The  final  section  considers 
the  topics  of  hidden  depths  and  the  variation  of  bottom  reflection  loss 
with  changes  in  sediment  type.  The  general  thrust  of  all  of  the  work 
reported  in  this  chapter  is  to  determine  the  importance  of  various  sub- 
bottom parameters  and  to  work  toward  a useful  description  of  the  sub- 
bottom which  requires  minimal  Information. 

Chapter  III  considers  the  sensitivity  of  propagation  loss  to  variations 
in  bottom  reflection  loss  using  a single  parameter  phenomenological  bottom 
loss.  Several  existing  propagation  loss  models  have  been  used  in  this 
study.  Major  questions  which  are  addressed  by  this  study  include  the 
delineation  of  regions  of  sensitivity  and  Insensitivity  to  bottom  reflec- 
tion loss,  as  well  as  the  way  this  separation  is  affected  by  soiorce  depth 
and  profile  type. 

Chapter  IV  considers  two  aspects  of  the  range  variable  bathymetry 
problem,  slope  enhancement  dijring  propagation  over  a continental  slope 
and  the  rough  interface  problem.  The  first  of  these  is  addressed  by 
considering  the  results  of  the  analysis  of  acoustic  data  taken  from  a 
recent  exercise  as  well  as  the  results  of  model  studies.  The  range  of 
bottom  irregularities  which  fall  into  the  category  of  bottom  roughness  is 
considered  as  well  as  various  approaches  to  accounting  for  the  associated 
acoustical  effects. 

Chapter  V considers  various  techniques  for  extracting  information 
about  subbottom  sediment  structure,  particularly  velocity  gradients,  from 
profiling  data.  This  chapter  constitutes  a particularly  concise  and  use- 
ful review  of  work  which  has  been  reported  and  is  not  well  known  in  the 
acoustics  conmmnity. 

Chapter  VI  serves  to  docvunent  the  work  carried  out  during  this 
reporting  period  on  the  development  of  normal  mode,  bottom  loss,  and  ray 
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trace  models  which  have  been  designed  with  a view  toward  the  solution  of 
bottom  interaction  problems. 


In  Chapters  II  through  VI,  figures,  tables,  equations,  and  references 
are  numbered  sequentially  within  each  chapter.  References  are  listed  at 
the  end  of  each  chapter. 

D.  Summary  of  the  Principal  Results 

This  section  contains  a brief  listing  of  the  major  results  obtained 
during  the  course  of  this  year's  research.  For  convenience,  these  results 
are  itemized  according  to  the  chapters  of  this  report  in  which  the  asso- 
ciated research  is  described. 

There  are  two  points  concerning  these  results  which  are  important  to 
appreciate:  (l)  all  conclusions  were  obtained  and  are  applicable  only 

over  some  range  of  conditions  and  under  circimistances  which  cannot  be 
detailed  in  this  section,  and  (2)  they  represent  only  the  results  of  the 
research  discussed  in  this  report  and  are  not  conclusions  of  this  total 
(ongoing)  research  program.  The  second  of  these  restrictions  is  one  that 
will  be  removed  in  time,  whereas  the  first  simply  reflects  the  fact  that 
results  taken  out  of  context  can  be  misleading. 

1.  Bottom  Loss  Studies 

The  presence  of  a basalt  basement  under  relatively  thin  clay 
sediments  containing  an  appreciable  sound  speed  gradient  can  lead  to 
unexpectedly  large  Increases  in  bottom  loss  over  a narrow  angular  range. 
This  loss  is  due  to  the  presence  of  Stoneley  waves. 

The  Stoneley  wave  loss  peak  occurs  at  the  low  grazing  angles 
which  are  crucial  in  determining  long  range  propagation. 

Attenuation  gradients  within  the  subbottom  can  be  important  in 
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determining  bottom  loss.  Ignoring  attenuation  gradients  can  lead  to 
errors  of  several  decibels  in  computed  bottom  loss. 


Density  gradients  within  the  subbottom  have  only  small  effects 
over  the  low  grazing  angle  regime  of  concern  to  long  range  propagation. 

The  effects  of  density  gradients  can  be  almost  wholly  accounted 
for  by  simply  using  the  local  density  values  at  each  interface  and  ignor- 
ing the  continuous  changes. 

The  sound  speed  gradient  within  the  subbottom  is  the  controlling 
parameter  in  determining  the  hidden  depth. 

The  hidden  depth  is  foimd  to  be  several  soimd  wavelengths  below 
the  ray  turning  depth.  This  result  is,  for  practical  purposes,  indepen- 
dent of  the  precise  definition  of  hidden  depth  so  long  as  it  is  defined 
in  terras  of  a fractional  change  in  bottom  loss. 

The  sensitivity  of  bottom  loss  to  changes  in  sediment  type  is 
best  studied  via  the  parameterization  of  density,  sound  speed,  and 
attenuation  by  porosity. 

This  parameterization,  the  results  of  which  have  been  used  in 
this  study,  has  been  carried  out  in  the  literature.  It  will  be  important 
to  extend  the  existing  parameterization,  developed  only  for  surficial 
sediments,  to  depth  in  the  sediment. 

2.  Propagation  Loss  Sensitivity  Studies 

Bottom  interaction  effects  on  propagation  loss  divide  the  water 
column  into  three  regions  based  on  source  and  source  conjugate  depths.  Propa- 
gation loss  has  a regular  dependence  on  critical  angle  in  the  central  region. 

Outside  this  region  of  regularity  (near  the  boundaries),  bottom 
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influence  Is  strongly  dependent  on  receiver  depth  and  is  more  sensitive 
to  bottom  loss  the  nearer  the  receiver  is  to  the  siirface. 

3 • Effects  of  a Range  Variable  Environment 

For  a shallow  source,  both  data  and  theoretical  calculations 
show  that  slope  enhancement  is  sensitive  to  receiver  depth,  at  least  for 
North  Atlantic,  double  minimum  type  profiles. 

When  slope  enhancement  is  potentially  present  geometrically,  it 
will  be  quite  sensitive  to  the  bottom  loss  near  the  upper  lip  of  the 
continental  slope. 

The  absence  of  enhancement  for  a deep  receiver  and  a shallow 
source  located  over  the  continental  shelf  is  associated  with  multipath 
conversion  process  which  severely  constrains  the  coupling  between  source 
and  receiver  for  this  geometry. 

The  definition  of  the  scope  of  the  rough  surface  problem, 
considering  the  length  and  height  scales  for  which  rox;ighness  will  be  impor- 
tant, as  well  as  probable  methods  for  obtaining  information  about  ocean 
bottom  roughness,  is  quite  sensitive  to  frequency  and  will  disappear  en- 
tirely in  the  very  low  frequency  (VLF)  region  submerging  litself  in  the 
(deterministic)  sloping  bottom  problem. 

The  most  promising  method  for  accoiintlng  for  roughness  effects 
on  propagation  is  the  smoothed  boundary  condition  approach. 

Rough  surface  scattering  theory  will  be  useful  primarily  for 
obtaining  estimates  for  the  amovint  of  scattering  loss  the  coherent  field 
will  suffer  per  bounce,  but  not  for  actually  carrying  out  quantitative 
propagation  calculations. 
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II.  BOTTOM  loss  STUDIES 


This  chapter  consists  of  four  Independent  sections  devoted  to  different 
aspects  of  the  problem  of  the  sensitivity  of  bottom  loss  to  variations 
of  subbottom  composition.  This  approach  is  part  of  a consistent  pattern 
of  attacking  the  bottom  interaction  problem  both  from  a propagation  view- 
point with  certain  classes  of  soiuid  speed  profile,  sensor  geometries,  etc., 
and  from  a more  local  viewpoint  using  bottom  reflection  loss  as  a measure 
of  the  effects  of  subbottom  variations. 

The  first  section  deals  with  the  effects  of  a solid  substrate  (basalt 
basement)  on  bottom  loss.  This  work,  the  beginnings  of  which  were  noted 
in  the  final  report  from  the  first  year's  work  in  this  program  (Hawker  et 
al.^),  deals  with  a loss  mechanism  which  can  significantly  affect  low  angle 
bottom  loss. 


The  second  and  third  sections  of  this  chapter  concern,  respectively, 
the  sensitivity  of  bottom  loss  to  the  variation  of  attenuation  and  density 
with  depth  in  the  sediment.  These  studies  are  directed  toward  establish- 
ing the  mlnimiam  amount  of  Information  concerning  subbottom  makeup  which  is 
required  for  accurate  prediction  of  either  propagation  loss  or  bottom  loss. 

The  fourth  section  of  this  chapter  concerns  the  problems  of  hidden 
depths  (the  depth  of  significant  acoustic  penetration)  and  sensitivity  to 
variations  in  sediment  type.  The  work  in  this  section  has  already  been 
reported  \inder  this  contract  (Hawker,  Focke,  and  Anderson  );  thus,  only 
a brief  synopsis  of  it  is  given  here. 

A.  The  Influence  of  Stoneley  Waves  on  Bottom  Reflection  Loss 


This  section  considers  the  effects  of  a particular  type  of  Interface, 
or  sxirface,  wave  known  as  a Stoneley  wave  on  bottom  reflection  loss.  These 
waves  can  exist  at  a solid-solid  bovindary  as  well  as  at  a fluid-solid 
boundary.  The  particular  Interface  of  concern  here  is  the  one  that  divides 
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sediments  (fluids)  from  the  basement  or  substrate  (solid),  usually  taken 
to  be  basalt. 

Typical  values  for  the  compressional  and  shear  wave  speeds  in  basalt 
are  5700  m/sec  and  2700  m/sec,  respectively.  Such  values  will  result  in 
compressional  and  shear  critical  angles  of  approximately  7^°  and  55° j 
respectively.  Given  the  constraints  on  grazing  angles  imposed  by  the  range 
of  propagation  geometries  of  interest,  it  is  clear  that  even  the  shear  wave 
critical  angle  will  not  constitute  an  important  effect  to  long  range 
propagation. 

At  first  sight,  then,  it  would  appear  that  the  relevant  effects  of 
the  basement  would  be  limited  to  a slight  Increase  in  bottom  loss  (rela- 
tive to  a perfect  reflector)  due  to  absorption  (nonperfect  reflection 
below  the  shear  wave  critical  angle).  This  is,  in  fact,  the  situation  for 
the  case  of  a solid  such  as  basalt  directly  underlying  the  water  column. 
However,  when  a sediment  layer,  either  a high  porosity  material  such  as 
clay  or  silt.  Intervenes  with  a positive  velocity  gradient  or  a low 
porosity  material  such  as  sand,  it  becomes  possible  to  excite  Stoneley 
waves  at  the  sediment -basalt  interface.  These  waves  (when  absorption  is 
present)  add  a new  loss  mechanism  which  can  result  in  large  increases  in 
bottom  loss  even  at  very  low  angles. 

The  effects  of  a solid  substrate  (basalt)  on  bottom  reflection  loss 
and  the  examination  of  the  various  wave  mechanisms  will  be  presented  in 
the  next  two  sections. 

All  bottom  loss  calculations  displayed  in  this  chapter  were  made 

using  the  ARL;UT  developed  bottom  loss  model  BOTLOSS,  described  in  detail 

3 

by  Hawker  and  Foreman. 

1.  The  Effects  of  a Solid  Substrate  on  Bottom  Reflection  Loss 


Figure  II-l  shows  schematically  the  bottom  model  to  be  considered. 
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PLANE 


WAVE  WAVE  (SV) 


Figure  n-i 

Schematic  Illustration  of  the  Wave  Structure  for  a 
Single  Fluid  Sediment  Layer  Overlying  a Solid  Substrate 
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AS- 77- 508 


The  water  is  a isovelocity  fluid  supporting  incoming  and  outgoing  plane 
waves,  the  amplitudes  of  which  are  related  by  the  complex  reflection 
coefficient  R.  The  sediment  can  have  a sound  speed  gradient  but  is 
otherwise  homogeneous.  The  substrate  half-space  is  permitted  to  support 
shear  waves  and  is  completely  homogeneous. 

The  bottom  loss,  -10  log,„  (RR*),  for  a clay  sediment  layer 

-1 

having  a sound  speed  gradient  of  1.2  sec  is  shown  in  Figs.  II-2  through 
II-^+  for  frequencies  of  25  Hz,  50  Hz,  and  100  Hz.  There  are  several 
characteristic  featiires  of  these  curves  which  should  be  noted. 

(1)  The  combination  of  a water  sound  speed  of  15^0  m/sec  and 
a substrate  compressional  wave  speed  of  5700  m/sec  produces  the  critical 
angle  effect  evident  at  approximately  7^°. 

(2)  There  is  a shear  wave  critical  angle  at  55°  resulting  from 
a substrate  shear  wave  speed  of  2700  m/sec. 

(3)  There  is  an  absence  of  the  shear  wave  critical  angle 
feature  at  55°  and  the  retention  of  the  compressional  wave  critical  angle 
effect  in  the  fluid  substrate  cases. 

(4)  There  is  a large  peak  in  loss  at  approximately  17°  in  all 
three  cases. 

(5)  The  absence  of  the  low  angle  j)eak  in  the  fluid  substrate 

cases. 

(6)  The  angular  location  of  the  low  angle  peak  is  nearly 
frequency  independent. 

Of  these  six  items,  (l)  through  (3)  are  well  understood  and 
present  no  difficulties  of  interpretation.  The  very  existence  of  a 
large  bottom  loss  peak  at  low  angles,  (4),  well  removed  from  the  shear 
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FIGURE  II-2 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
FOR  A 100  m CLAY  LAYER  AT  25  Hz 
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FIGURE  IE-4 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
FOR  A 100  m CLAY  LAYER  AT  100  Hi 


wave  critical  angle  and  clearly  related  to  the  solid  aspects  of  the  substrate 
(5),  not  to  mention  the  frequency  independence,  (6),  is  evidence  of  some 
additional  and  different  loss  mechanism. 


Figures  II-5  and  II-6  show  additional  evidence  for  such  a 
mechanism.  In  these  figures  the  bottom  loss  versus  grazing  angle  is 
depicted  as  a function  of  the  thickness  of  the  clay  layer.  The  large 
"ridge"  composed  of  the  anomalous  low  angle  loss,  absent  in  the  fluid 
substrate  case,  is  seen  to  be  a strong  function  of  the  layer  thickness. 

The  minor  ridge  seen  in  Fig.  II-6  can  be  shown  to  be  located  at  the  depth 
at  which  the  turning  depth  equals  the  layer  thickness.  A comparison  of 
Figs.  II-5  and  II-6  shows  that  the  anomalous  low  angle  loss  peaks  occur  at 
depths  just  below  the  depths  defined  by  this  turning  depth  related  ridge. 

Thus  far  we  have  considered  only  high  porosity  (clay)  sediments 

with  a sound  speed  gradient.  Figure  II-T  shows  bottom  loss  versus  grazing 

angle  for  a sand  layer  50  m thick.  The  sand  has  been  taken  to  have  zero 

atten\iation  simply  to  display  vividly  the  critical  angle,  at  24°  as 

well  as  the  substrate  shear  wave  critical  angle  9 at  57°.  The  large 

s 

peak  in  loss  at  20°,  below  the  critical  angle,  corresponds  to  the  low 
angle  peaks  seen  previously  for  clay  sediments.  In  the  sand  layer  case, 
however,  the  sound  speed  gradient  has  been  taken  to  be  zero. 

Figure  II-8  shows  the  same  situation  as  Fig.  II-T  except  that  the 
layer  is  75  id  thick.  Unlike  the  clay  layer  case  for  which  the  peak  shifts 
location  with  layer  thickness  but  maintains  a large  amplitude,  the  iso- 
velocity sand  layer  case  shows  no  location  change  but  a reduced  amplitude. 

A series  of  calculations  made  with  a fixed  thickness  clay  layer  shows  that 
the  low  angle  peak  also  shifts  location  with  sound  speed  gradient  and 
disappears  in  the  zero  gradient  case.  These  various  aspects  of  the  low 
angle  peaks  will  be  fully  explained  by  the  Stoneley  wave  mechanism  to  be 
introduced  in  the  second  subsection. 

Finally,  we  consider  a sequence  of  sediment  types  having  a 
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FREQUENCY,  50  Hi 
LAYER  MATERIAL:  CLAY 
SOUND  SPEED  GRADIENT:  1.0  see"’ 
SUBSTRATE:  SOLID 


figure  IE-5 

BOTTOM  LOSS  versus  GRAZING  ANGLE  arl  - ut 

versus  LAYER  THICKNESS  FOR  A CLAY  LAYER  kIhVrVo 

OVERLYING  A SOLID  SUBSTRATE  6 - 28  - 76 
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BOTTOM  LOSS  - dB 


FREQUENCY:  50  Hz 
LAYER  MATERIAL:  CLAY 
SOUND  SPEED  GRADIENT:  1.0  sec"’ 
SUBSTRATE:  PERFECT  REFLECTOR 


FIGURE  II-6 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
versus  LAYER  THICKNESS  FOR  A CLAY  LAYER 
OVERLYING  A PERFECT  REFLECTOR 
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FREQUENCY;  50  Hz 

SEDIMENT:  SAND  WITH  ZERO  ATTENUATION 
SUBSTRATE:  BASALT 
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BOTTOM  LOSS  versos  GRAZING  ANGLE  FOR  AN 
IDEAL  75  m ISOVELOCITY  SAND  LAYER  AT  50  Hi 


sound  speed  gradient  of  1.2  sec  ^ and  a thickness  of  100  m.  The  density, 

surficial  sound  speed,  and  attenuation  are  related  to  porosity  by  the 

T 

regression  equations  of  Hamilton  and  the  resulting  parameters  used  to 

compute  the  bottom  loss.  The  results  are  shown  in  Fig.  11-9^  which  is 

2 

taken  from  Hawker,  Focke,  and  Anderson. 

The  curving  ridge  of  bottom  loss  peaks  between  15°  and  25°  shows 
the  effect  on  the  low  angle  loss  peak  of  varying  the  sediment  type.  As 
would  be  expected  on  the  basis  of  Figs.  II-7  and  II-8,  the  sand  sediments 
(porosity  = to  do  not  show  an  anomalous  peak  for  a layer  of  this 

thickness.  The  peaks  are  most  pronounced  for  the  high  porosity  sediments, 
silts,  and  clays. 

2 . Stoneley  Waves  and  Bottom  Reflection  Loss 

Before  entering  into  a discussion  of  the  nature  of  Stoneley 
waves  at  a fluid-solid  interface,  it  will  be  useful  to  carry  forward 
the  heuristic  studies  of  the  previous  subsection  in  order  to  fiorther 
understand  the  characteristics  of  the  loss  mechanism. 

Figure  II-IO  shows  the  same  situation  as  in  Fig.  II-5  except  that, 
in  the  present  case,  the  sediment  (clay)  attenuation  as  well  as  the  sub- 
strate compressional  attenuation  have  been  set  to  zero.  The  17°  bottom 
loss  peak  is  still  evident  in  Fig.  II-IO;  however,  below  the  neighborhood 
of  this  peak,  there  is  perfect  reflection  (zero  loss).  Although  examples 
are  not  shown  here,  essentially  the  same  thing  occurs  if  any  two  of  the 
three  relevant  attenuations  are  set  to  zero  with  the  remaining  attenuation 
(the  substrate  shear  attenuation  in  the  present  case)  nonzero. 

Figure  II-ll  shows  the  result  of  setting  all  three  attenuations 
to  zero.  As  expected,  there  is  no  loss  below  the  shear  critical  angle 
(55°)  whereupon  the  loss  increases  and  remains  nonzero  up  to  90° • The 
results  shown  in  Figs.  II-IO  and  II-ll  show  that  the  mechanism  responsible 
for  the  low  angle  peaks  is  dependent  on  at  least  one  medium  composing 
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FIGURE  H-  11 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
WITH  ALL  ATTENUATIONS  ZERO 


the  interface  being  lossy.  These  peaks  are  not  simply  due  to  energy 
being  propagated  away,  which  is  the  mechanism  responsible  for  the  loss 
shown  in  Fig.  II-ll. 


Figures  11-12  and  11-13  show  the  phase  of  the  reflection 
coefficient  for  the  same  situation  used  in  Figs.  II -3,  II-IO,  and  II-ll; 
for  a solid  substrate.  Fig.  11-12;  and  a fluid  substrate.  Fig.  11-13.  It 
will  be  observed  that  in  the  neighborhood  of  the  17°  loss  peak  (denoted  by 
a cross  on  Fig.  11-12)  the  phase  is  rapidly  changing  and  goes  thro\;gh  an 
entire  extra  cycle  of  variation  relative  to  the  fluid  case. 

Figure  II-l4  shows  the  pressxire  amplitude  versus  depth  within 
the  same  100  m clay  layer  for  a plane  wave  incident  at  16.5°,  for  both 
fluid  and  solid  substrates.  The  fluid  case  displays  the  expected  behavior, 
oscillatory  above  the  turning  depth,  Z^,  and  quasiexponentially  decaying 
below  Z^.  The  solid  substrate  case  is  similar  except  the  field  shows  a 
growing  amplitude  with  depth  very  near  the  surface.  Figure  II- I5  is  the 
same  situation  as  Fig.  II-l4  except  that  the  grazing  angle  is  17 .5° > flie 
angle  at  which  the  peak  in  loss  occurs.  The  fluid  substrate  case  is 
essentially  unchanged  (the  turning  depth  is  slightly  deeper);  however, 
the  solid  case  is  substantially  different.  The  field  now  displays  a 
sharp  Increase  in  amplitude  as  the  interface  is  approached,  with  a near 
constant  behavior  above  Z^. 

Of  course,  since  the  substrate  is  taken  to  be  homogeneous,  the 
shear  and  compressional  fields  are  purely  linear  combinations  of  exponential 
solutions.  For  angles  below  the  shear  and  compressional  critical  angles, 
the  shear  and  compressional  wave  potential  are  respectively 

i(k  xcos0-cut) 

%=  A e ^ e , 

s 


l(k  xcosS-oit) 

$ = A e 
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FIGURE  n-I3 

REFLECTION  COEFFICIENT  PHASE  versus  GRAZING  ANGLE 
FOR  A FLUID  SUBSTRATE 


■'T'll 


- HidaO  6-7-76 

30 


6-7-76 


- Hidaa 


31 


2 2 2 2 

K = k cos  0 - k 
so  s 

2 ,2  2 „ ,2 

K = k cos  0 - k 
P o p 


k is  the  water  wave  number,  and  k and  k are  the  substrate  shear  and 
o ' s p 

compressional  wave  numbers.  The  amplitudes  and  are  constants. 

The  fields  therefore  decay  away  from  the  Interface  in  both 

directions  when  the  grazing  angle  is  near  the  angle  at  which  the  low 

angle  loss  peak  occurs.  This  behavior  is  characteristic  of  an  interface 

(or  surface)  wave  and  without  further  evidence  can  be  taken  to  prove 

that  a Stoneley  wave  is  present.  This  point  is  discussed  in  detail  by 
4 

Tolstoy  and  Clay  for  two  joined  Isovelocity  half-spaces  which  lead  to  a 
dispersion  equation  for  the  \inlque  phase  velocity  of  Stoneley  waves. 

If  we  ignore  for  the  moment  the  complexity  introduced  by  the 

layered  inhomogeneous  config\iration  being  treated  here,  we  can  understand 

how  there  must  be  a unique  angle  associated  with  the  excitation  of  Stoneley 

leaves.  Since  these  waves  can  have  only  one  phase  velocity,  v^^  and  this 

is  in  the  direction,  x,  parallel  to  the  interface,  they  can  exist  only 

when  the  plane  waves  have  an  equal  x component.  That  is,  c^/cos  9 must 

equal  v , thus  defining  an  angle  0. 
s 

A sufficiently  accurate  model  which  contains  the  basic  physical 
processes  can  be  obtained  by  considering  a fluid  half-space  overlying  a 
solid  half-space.  The  fluid  Is  taken  to  have  a sound  speed  which  decreases 
away  from  the  interface,  while  the  solid  is  assumed  to  be  isovelocity. 

The  neglect  of  the  water- sediment  interface  bo\inding  the  fluid  layer  is 
based,  in  part,  on  the  fact  thet  Fig.  11-15  shows  the  field  decaying 
quaslexponentially  away  from  the  sediment- substrate  interface  for  several 
sound  wavelengths.  The  nature  of  the  field  in  the  region  near  this  lower 
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interface  will  thus  be  effectively  isolated  from  the  presence  of  the  upper 
interface  which  may  therefore  be  removed. 


The  z-axis  is  taken  to  be  increasing  positively  downward  (into 
the  solid)  with  the  fluid-solid  interface  located  at  z=0.  With  the  fluid 
velocity  potential  denoted  by  and  the  compressional  and  shear  potentials 
in  the  solid  denoted  by  <I>g  and  we  have  the  usual  fluid-solid  relations, 


k 


b\|f^  dijfg 


^ ^ 


0=2' 


dxdz  ^ 


dJ 


dz^ 


(see,  for  example,  Brekhovsklkh) . 


The  fields  are  now  assumed  to  be  of  the  forms 


= A^f(z) 

/^p%i(ax-cot) 


222/2  222  2 

where  ic  -a  -kd  /c  =0,  k -a  -hu  /c  =0  and  the  amplitudes  A , A , B are 
p P s ' s ^ 1 2’  2 

unknown.  The  depth  function  f(z)  in  the  fliv'd  (O^^--®)  is  determined  by 

the  sound  speed  profile  c^(z)  and  a radiation  condition. 

If  a combination  of  amplitudes  can  be  found  for  some  horizontal 
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wavenumber  a such  that  the  continuity  conditions  can  be  satisfied,  then 
solutions  of  the  type  proposed  here  can  exist.  Substitution  of  the 
assumed  forms  for  and  \|(  into  the  continuity  conditions  clearly 

leads  to  three  simultaneous  homogeneous  linear  algebraic  equations  for 
Af,  A^,  and  These  equations  have  a nontrivial  solution  only  if 

their  determinant  vanishes.  Upon  requiring  this,  a single  (dispersion) 
equation  is  obtained  for  the  horizontal  wave  number  a,  or  the  associated 
phase  velocity  v 


v^  + ^ 

c ^ ^2  c f'(0) 

P s 


This  equation  can  be  rewritten  in  a slightly  more  perspicacious 
form  by  introducing  the  fluid  sound  speed  at  the  interface  c (z=0)=c 

I 2 ^ ^ 

and  the  function  e(v)  = u)f(0)/(vf  ’ (O) ).  /l-v  / 2 , to  obtain 

'^1 


In  the  case  e(v)=l,  this  equation  becomes  Identical  with  the  classical 
result  for  the  phase  velocity  of  Stoneley  waves  at  the  interface  between 

4 

two  isoveJoclty  media,  as  given  by  Tolstoy  and  Clay. 


Before  discussing  solutions  to  this  equation,  it  is  worth 
recapitulating  the  situation  thus  far.  Solutions  of  a particular  type, 
defined  by  the  forms  assumed  for  0^,  and  can  exist  in  the  two 
media  if,  and  only  if,  the  phase  velocity,  v,  obeys  the  dispersion 
equation  above. 

If  the  function,  e(v).  Itself  contains  an  arbitrary  constant, 
such  as  that  relating  incoming  and  outgoing  plane  waves  when  c^(z)  is 
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constant,  then  a solution  for  this  equation  is  obviously  possible  for 
all  values  of  v.  Thus,  for  the  plane  wave  reflection  case,  no  re- 
striction is  imposed  by  the  dispersion  equation  and  waves  can  exist  for 
all  phase  velocities  (wave  numbers).  However,  when  a single  solution, 
f(z),  is  sought,  no  arbitrariness  remains  and  only  discrete  values  of 
V are  possible . 

Since  we  have  already  assumed  that  only  the  sediment  properties 
near  the  fluid-solid  interface  are  of  importance,  it  is  perfectly  consis- 
tent to  approximate  the  linear  sound  speed  profile  used  in  the  bottom 
loss  calculations  of  subsection  1 by  the  pseudolinear  profile  Cj^(2)~c^(0)/ 
(l-Pz).  This  form  leads  to  the  differential  equation, 

f"  + (k^^-a^-k^^Pz)f  = 0 , 

for  f in  the  region  0^.  The  value  of  ^ is  taken  to  be  2g/c^(0)  where 

g is  the  desired  sound  speed  gradient  at  the  interface.  The  solutions 

2 2 

to  this  equation  are  the  Airy  functions  Ai(u)  and  Bi(u)  with  u={a  -k  ) 

2 -*2/3  2 1/3 

(k^  P)  ''^+(k^  P)  '""^z.  Since  we  wish  to  investigate  the  possibility  of 
solutions  which  decrease  in  amplitude  away  from  the  interface,  we  choose 
the  solution  f=Bi{u). 

The  resulting  equation  has  been  solved  numerically  for  the  phase 
velocity  V given  the  parameters  {g,  c , c , c , cu,  p /p_}.  Since  in  the 

X p S J.  ^ 

plane  wave,  layered  reflection  problem  that  is  being  approximated  here, 
the  X dependence  was  given  by  exi(ik^xcose),  we  must  have  Q!=k^  cos  9 = 
cn/v  = cu  cos  P/c^.  The  angle  6 thus  defined  will  be  the  angle  at  which  a 
boundary  wave  can  exist. 

The  results  of  such  numerical  calculations  are  given  in  Figs. 
II-I6  throxigh  II-I8,  which  show  the  predicted  Stoneley  wave  excitation 
angle  (solid  line)  and  the  angles  at  which  direct  bottom  loss  calculations 
show  the  low  angle  anomalous  peak  (circles).  Figure  II- I6  shows  these  two 
quantities  plotted  versus  layer  thickness  for  a clay  layer  having  a sotind 
speed  gradient  of  1.2  sec  , whereas  Fig.  II-IT  employs  the  sound  speed 
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FIGURE  H-  18 

LOCATION  OF  THE  LOW  ANGLE  BOTTOM 
LOSS  PEAK  versus  FREQUENCY 


gradient  (with  fixed  layer  thickness)  as  the  independent  variable. 
Agreement  is  found  to  be  excellent  in  both  cases. 

The  slight  frequency  dependence  shown  in  Fig.  II-I8  shows  that 
the  very  small  shift  of  the  loss  peak  with  frequency  shown  in  Figs.  II-2 
through  II-4  is  also  accounted  for  with  this  theory. 

These  calculations  provide  clear  evidence  that  the  low  angle 
bottom  loss  peak  discussed  in  the  first  subsection  of  this  chapter  is 
indeed  due  to  Stoneley  waves.  Several  additional  qualitative  features 
are  also  accounted  for  by  this  mechanism.  For  instance,  the  requirement 
for  at  least  one  nonvanishing  absorption,  as  evidenced  by  Figs.  II-5, 
II-IO,  and  II-ll,  is  now  understood.  Since  the  Stoneley  wave  propagates 
along  the  interface,  it  can  carry  away  no  energy  when  it  is  excited  by  a 
plane  wave.  However,  when  absorption  is  present,  dissipation  occurs 
along  the  propagation  path  and  loss  results.  Other  features  such  as  the 
density  dependence  of  the  location  of  the  peak  can  also  be  explained  by 
the  Stoneley  wave  mechanism. 

5 . Practical  Implications 

The  basic  mechanism  causing  the  anomalous  low  angle  bottom  loss 
peaks  has  been  explained  and  understood.  It  is  now  meaningful  to  turn 
our  attention  to  the  question  of  the  practical  implications  of  this  loss 
mechanism. 

There  are  several  factors  evident  in  the  work  reported  here 
which  tend  to  increase  the  potential  importance  of  this  loss  mechanism  to 
propagation  problems: 

(1)  the  presence  of  a large  loss  peak  in  the  low  angle  regime, 

0°  to  25”, 

(2)  the  possibility  of  exciting  Stoneley  waves  even  for  fairly 
thick  sediment  overburdens,  especially  in  case  of  clay  sediments,  and 
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(5 ) the  apparent  broadening  of  the  peak  width  with  decreasing 

frequency. 


There  are,  on  the  other  hand,  factors  whose  impact  on  the  Stoneley 
wave  loss  peaks  are  at  present  unknown  and  are  the  subject  of  continuing 
investigation: 

(1)  the  influence  of  shear  waves  in  the  sediments  themselves, 

(2)  the  effects  of  lateral  variations  in  sediment  properties 
and  layer  thickness,  and 

(5)  the  effects  of  Interface  irregularities  and  density 
inhomogenities . 

In  an  ox)en  ocean  measurement  of  bottom  loss,  data  scatter  is 
frequently  so  great  as  to  preclude  direct  determination  of  the  presence 
or  absence  of  relatively  narrow  bottom  loss  peaks.  Nevertheless,  such 
peaks  could  manifest  themselves  in  the  data.  In  particular,  a large 
though  comparatively  narrow  loss  peak  at  a low  grazing  angle,  not  directly 
resolvable,  could  result  in  a general  increase  in  measured  bottom  loss. 

It  is  plausible  that  measurements  showing  high  bottom  loss  (>5  dB)  at  low 
angles  (<10®)  for  clay  sediments  are  influenced  by  such  loss  peaks. 

B.  Sensitivity  of  Bottom  Loss  to  the  Depth  Variation  of  Attenuation 

In  this  section,  the  effects  of  various  treatments  of  attenuation 
on  bottom  loss  curves  are  demonstrated.  The  program  BOTIOSS,  described 
in  Chapter  V,  was  used  to  compute  all  the  bottom  loss  curves  in  this 
chapter. 

The  shape  of  the  bottom  loss  curve  is  determined  primarily  by  the 
so\ind  speed  profile.  Attenuation  increases  bottom  loss  for  angles  where 
compressional  waves  enter  the  bottom.  An  attenuation  gradient  Increases 
the  bottom  loss  even  more.  This  effect  will  be  examined  for  a clay 
bottom.  Some  preliminary  judgments  concerning  the  validity  of  current 
attenuation  treatments  will  be  made  using  the  best  current  attenuation 
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data  for  comparison. 


In  situ  measurements  of  attenuation  versus  depth  have  recently  been 

5 

reviewed  by  Hamilton,  who  discusses  basic  sand  and  clay  bottom  types  and 
variations  in  the  behavior  of  atteniiatlon  with  depth.  Figure  II-I9 
illustrates  the  conclusions  of  his  review.  Hamilton's  collected  data  for 
clay  bottoms  cover  primarily  the  upper  200  ra.  At  depths  greater  than  400 
all  sediments  seem  to  have  the  same  attenuation  profile.  Shirley^  has 
suggested  the  form  of  the  silt-clay  attenuation  profile  between  200 
and  UOO  m. 

Clay  bottoms  are  characterized  by  low  reflectivities  and  moderate 
velocity  gradients  (l.O  to  I.5  sec  ^) . Attenuation  increases  with  depth 
to  between  200  and  400  m.  The  velocity  at  the  sediment-water  interface 
is  low.  with 


0.90  < 


^ 1.05 


Table  II-l  presents  the  input  parameters  used  for  the  bottom  loss  curves 
in  this  chapter. 

The  results  for  an  isovelocity,  Isodenslty  layer  300  m thick  are 
shown  in  Fig.  11-20.  With  no  attenuation,  no  bottom  loss  occurs  until  the 
substrate  shear  critical  angle  at  62°  is  reached.  A constant  attenuation 
leads  to  nonzero  bottom  loss  below  62°  and  a small  increase  in  loss  above 
62°.  Attenuation  linearly  increasing  with  depth  causes  an  additional  in- 
crease in  bottom  loss  beyond  the  constant  attenuation  core.  The  attenua- 
tion gradient  exaggerates  the  high  loss  areas  of  the  curve. 

Figure  11-21  shows  the  results  for  a 300  m layer  with  linearly 
Increasing  velocity  and  density.  Note  that  the  change  in  sound  speed 
treatment  from  constant  to  linearly  increasing  greatly  alters  the  shape 
of  the  curve.  The  comments  for  Fig.  11-20  concerning  attenuation  effects 
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also  apply  to  Fig.  11-21.  It  is  unfortunate  that  sound  speed  profiles 
below  the  top  few  meters  of  sediment  are  very  often  poorly  known. 

Figure  11-22  demonstrates  the  effect  of  various  treatments  of 
attenuation.  Sound  speed  and  density  are  linearly  increasing  with  depth. 
The  parameters  listed  in  Table  II-l  are  used  for  a layer  depth  of  500  ra. 
The  modified  exponential  profile  and  sound-speed-linked  profile  are 
described  in  the  last  section  of  this  chapter.  All  of  the  attenuation 
profiles  used  are  sketched  in  Fig.  11-25.  The  pseudolinear  profile  con- 
tained a singularity  at  a shallow  depth  and  could  not  be  used.  Hamilton's 
clay  data  region  is  indicated  for  comparison  purposes. 

The  treatments  fall  into  two  groups.  The  constant  and  sound-speed- 
linked  attenuation  profiles  give  very  much  the  same  results.  In  the 
sound- speed- linked  case,  attenuation  is  proportional  to  c^/c^*  This 
causes  attenuation  to  decrease  with  depth  and  makes  it  inconsistent  with 
the  observed  behavior. 

The  linear  and  modified  exponential  profiles  give  similar  bottom  loss 
curves.  The  linear  profile  bottom  loss  curve  is  1 to  5 dB  higher  than 
the  modified  exponential  profile  curve,  mainly  in  the  higher  loss  regions 
of  the  curve.  Either  of  these  treatments  is  consistent  with  the  current 
attenuation  data. 

Since  attenuation  is  cumulative  over  the  path  of  the  compressional 
wave  in  the  sediment  and  the  penetration  of  waves  into  the  bottom  is  de- 
termined by  their  incidence  angle,  the  use  of  an  averaged  attenuation 
over  the  entire  layer  gives  a different  curve  than  the  use  of  a linear 
or  other  attenuation  profile  for  thick  layers. 

In  Fig.  11-24,  note  the  higher  bottom  loss  for  angles  letc,  than  the 
critical  angle  at  52°.  Around  10°  the  difference  is  1.5  dB,  with  the 
averaged  attenuation  curve  higher  than  the  linear  attenuation  profile 
bottom  loss  curve.  After  the  52°  sediment  critical  angle,  the  linear  and 
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FIGURE  ir-24 
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averaged  atteniiation  profiles  are  nearly  Identical.  The  sound  speed  and 
density  values  and  gradients  are  as  in  Table  II-l. 

Further  attenuation  studies  now  in  progress  include  a study  of  sand 
bottoms,  a study  of  the  effects  that  the  change  of  attenuation  gradient 
from  positive  to  negative  have  on  bottom  loss  in  the  clay  case,  and  studies 
of  attenuation  effects  in  more  realistic  multilayered  bottoms  of  various 
types.  The  results  of  these  studies  vill  be  discussed  in  a future  report. 


C.  Sensitivity  of  Bottom  Loss  to  Subbottom  Density  Gradients 

The  purpose  of  this  section  is  to  report  the  results  of  a sensitivity 
study  on  the  effects  of  sediment  density  gradients  on  bottom  loss.  This 
section  is  divided  into  four  subsections.  The  first  subsection  will  dis- 
cuss theoretical  treatments  of  the  problem  including  the  results  obtained 
from  bottom  loss  models  to  support  the  conclusions  of  the  theoretical 
treatments.  The  second  and  third  subsections  will  classifjf  the  magnitude 
of  the  density  gradient  effects  as  to  sediment  type  and  grazing  angle 
regimes.  The  foiarth  s\ibsection  summarizes  the  results  and  conclusions 
of  this  study. 

In  this  sensitivity  study  it  was  necessary  to  adopt  two  models,  one  being 
a bottom  loss  model  and  the  other  being  a sediment  structure  (geoacoustic) 
model.  The  model  used  to  evaluate  the  bottom  loss  as  a function  of  grazing 
angle  was  BOTLOSS,  a plane  wave  reflection  coefficient  model  based  on 
numerical  integration.  BOTLOSS  computes  the  bottom  loss  as  a ihinction 
of  grazing  angle  for  a multilayer  fluid  sediment  structure  overlying  a 
solid  substrate  by  numerically  Integrating  the  wave  equation  through 
the  different  layers  and  applying  appropriate  boundary  conditions  at 
the  layer  interfaces. 

1.  Background 

The  wave  equation  in  the  jth  subbottora  layer  when  a density 
gradient  in  the  vertical  direction  is  present  1s 


J J J 


1 a 

— p P 

Cl/  \ J 


c . ( z ) dt 
J 


where 


If  Pj(xzt)  is  expressed  as 


. , ik  COS0X 

, , - lent  o / \ 

P.  (x,z,t)  = e e P,U) 

J J 


the  following  equation  for  p.(z)  is  obtained 

J 


‘^^Pl  r ? 2 2 T 

^ (z)  + r.(z)— ^ + [k.  (z)  - k cos  e]p.(z)  = 0 

dz  J az  J o J 


(1) 


The  angle  9 is  the  grazing  angle  at  the  boundary  between  the  water  and 
subbottom  layers  and  k^  is  the  wave  number  in  the  water. 

A useful  form  of  Eq.  (l)  is  obtained  by  reducing  it  to  normal  form. 

This  is  accomplished  by  writing  p.(z)  as 

J 

z 

Pj(z)  = Xj(z)  exp|-  r(z')dz'|  , 


where  z^  is  an  arbitrary  reference  depth. 


If  this  equation  is  used  in  Eq.  (l),  the  following  differential 
equation  for  Xj(z)  is  obtained. 


2 


(2) 
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with 


P.(z)  = 

c) 


p.(z) 


1/2 


X.(z) 

J 


(5) 


An  expression  for  the  reflection  coefficient  in  terms  of  the  wave 
equation  in  the  first  subbottom  layer,  p^(z),  may  be  obtained  by  an 
application  of  the  continuity  conditions  at  the  water- sediment  interface. 
This  procedure  yields 


-(po/Pl)^P'A^o 

“ +(Po/Pi7p'i/Pi^o 


(M 


where  K = k sin  9,  p is  the  water  density,  and  p is  the  density  of  the 
o o o ./  / 

first  subbottom  layer  at  the  water- sediment  interface,  i.e.,  z=0.  Equa- 
tions (2)  through  (4)  will  be  used  in  the  upcoming  theoretical  discussion. 
For  details  concerning  the  actual  numerical  solution  of  the  wave  equations 
and  the  boundary  conditions,  see  Hawker  and  Foreman.^ 


The  geoacoustic  model  used  for  the  bottom  in  this  sensitivity  study 
is  a simple  one.  The  subbottom  structure  is  modeled  as  a relatively  thin 
fluid  sediment  layer  (100  to  JOO  m)  having  a positive  velocity  gradient, 
overlying  an  Isovelocity  solid  substrate  half-space,  in  which  the  substrate 
is  treated  as  a true  solid  supporting  shear  waves  as  well  as  compressional 
waves  (see  Fig.  11-25).  This  layering  structure  models  the  abyssal  plain 
environment  reasonably  well. 


The  examples  used  in  this  study  will  deal  with  three  different 
sediment  types:  clay,  silt,  and  sand.  The  acoustic  parameters  of  the 
sediments  that  were  used  in  BOTLOSS  are  listed  in  Table  II-2. 


A feeling  for  the  effects  of  density  gradients  may  be  obtained  by 
an  examination  of  Eqs.  (2)  and  (5). 


FIGURE  11-25 

LAYER  CONFIGURATION  AND  VELOCITY  STRUCTURE 


ARL  - UT 
AS-76- 1213 
SRR  - DR 


TABLE  I 1-2 


SEDIMENT  ACOUSTIC  PARAMETERS 


Water  Density  = 1-055 


Water  Sound  Velocity  = 15^  m/sec 


Clay 

Silt 

Sand 

Velocity  at  top 
of  layer 

1526.6  m/sec 

1518.3  m/sec 

1723 .8  m/ sec 

Density  at  top 
of  layer 

1.27  g/cm^ 

1.68  g/cm^ 

2.08  g/cm^ 

Attenuation 

0.0029  dB/m 

0.012  dB/m 

0.024  dB/m 

Substrate  Density  =2.6  g/cm^ 
Compressional  Velocity  = 5700  m/sec 
Compressional  Attenuation  = 0.0015  dB/m 
Shear  Velocity  = 2700  m/sec 
Shear  Attenuation  = 0.01  dB/m 
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dz 


2 <j(z)  + 


k.^(z)  - k ^ cos  ^0  K 
J o 


(z)  = 0 


k.^(z)  = k.2(z)  . ir^.(z)  - |r- .(z) 


p.(z)  = 

J 


p.(z) 

— U-, 


p^(z  ) 

■ J o . 


1/2 


Xj(z) 


Equation  (2)  reveals  that  one  effect  of  a density  gradient  is  to  alter 

the  sound  speed  profile  from  c.(z)  to  c.(z),  where 

<3  J 


c (z)  = O)  [k^(z)  - Jr  ^(z)  - Jr  '(z)]’^^^ 


2‘  j 


Thus  a positive  density  gradient  causes  c.(z)  to  increase  relative  to 

J 

the  case  p'  =0.  Hence,  a positive  density  gradient  will  introduce  addi- 
J 

tlonal  refractive  effects  beyond  those  caused  by  a positive  sound  speed 
gradient.  This  modification  to  the  sound  speed  profile,  however,  is 
very  small  for  the  values  of  p'/p  expected  to  be  encountered  in  the  ocean. 


A very  simple  estimate  of  the  effect  of  a density  gradient  on  the 
sound  speed  profile  may  be  made  on  the  basis  of  a model  in  which  c (z)  is 

Y Z ^3 

assumed  constant  and  pj(z)=p(0)  e with  Y^O.  In  this  case 

ri(z)=-Y  and  r '(z)=0.  Hence,  c .(z)=<ni  k ^-Y^/^]  =C.[l-Y^/4k,^J 

j J J J J 

Typical  values  of  p'/p  deep  sea  sediments  would  be  on  the  order  of 

10  ^ m ^ (see  Hamilton  . Hence, Y^/4  k ^=(Yc/4'irf)^<1.5><10  ^/f^. 

O2  Jg 

Therefore,  for  frequencies  above  10  Hz,  Y /4k  «1;  this  implies  that 

J 

such  a modification  would  be  negligible. 


A higher  order  density  gradient  effect  is  seen  in  Eq.  II-J.  Equation 
II-5  indicates  that  an  added  modification  Introduced  by  the  density  gradient 
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l/2 

is  through  the  multiplicative  factor  [p(z)/p(z^)]  ' , This  multiplicative 

factor  represents  a modification  to  the  pressure  field  that  Increases 
with  depth  in  the  presence  of  a positive  density  gradient.  This  factor, 
although  usually  not  too  different  from  1,  does  represent  an  effect  of 
higher  order  than  the  modification  introduced  into  the  velocity  profile. 

The  two  modifications  just  discussed  are  due  to  added  terms  in  the 
wave  equation  which  are  present  when  density  gradients  are  considered. 
Another  "hidden  effect”  arises  from  the  boundary  conditions  that  the 
pressure  wave  function  must  satisfy  across  layer  boundaries.  When  the 
density  is  variable,  it  is  only  the  densities  at  the  interfaces  which 
are  used  in  the  continuity  conditions.  Since  the  interface  reflection 
and  transmission  coefficients  are  dependent  upon  the  density  ratios  across 
the  Interface,  it  does  not  seem  unlikely  that  this  "hidden  effect"  can 
be  dominant  in  many  circumstances.  Later  in  this  section  some  examples 
will  be  given  showing  that  the  effect  of  a density  gradient  can  be 
handled  quite  adequately  by  using  Interface  densities  in  the  boundary 
conditions . 

Since  the  models  used  in  this  investigation  explicitly  assume 
positive  sound  speed  gradients,  it  is  convenient  to  consider  low  grazing 
angle  and  high  grazing  angle  density  effects  separately.  In  the  lew- 
angle  region  where  refractive  effects  keep  the  incident  energy  away 
from  the  substrate,  the  density  gradient  effects  will  be  mainly  due  to  the 
presence  of  the  density  gradient  in  the  wave  equation.  In  the  high  angle 
regions  where  significant  energy  encomters  the  bottom,  it  is  expected 
that  the  hidden  "boundary  condition  effects"  will  become  dominant. 

2 . Low  Grazing  Angle  Effects 

To  derive  an  expression  which  is  useful  in  the  low  angle  region 
and  Indicates  the  effect  of  a density  gradient  on  the  magnitude  of  the 
reflection  coefficient,  it  is  necessary  to  adopt  a tractable  mathematical 
model.  The  following  development  will  consider  a two  half-space  problem. 
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The  first  half-space  will  be  isovelocity  and  the  second  will  be  explicitly 
assumed  to  have  a positive  velocity  gradient  and  an  exponential  density 
variation  given  by  p(z)=p(z^)  e^^.  The  reason  that  this  mode]  is  ex- 
pected to  work  reasonably  well  for  the  c\arrent  sediment  model  is  that, 
at  low  angles,  rays  incident  on  the  bottom  have  shallow  turning  depths 
due  to  velocity  gradients.  These  turning  depths  are  usually  shallow 
enough  such  that  the  exponential  density  variation  is  essentially  linear 
and  the  substrate  at  deeper  depths  is  effectively  absent  from  the  problem. 

Consider  two  half-spaces,  the  upper  being  isovelocity  with  a 
sound  speed  and  density  of  c^  and  and  the  lower  having  an  arbitrary 
positive  velocity  gradient  and  an  exponential  density  variation  [p(z)=p 
(z^)  e^^].  This  being  the  case,  the  pressure  variation  with  depth  in 
the  lower  layer  is  given  by 


3^(z)  = • 


P3^(z)' 


1/2  Yz/2 

X^(z)=e  x^(z) 


where  X2^(z)  satisfies 

£ 

■ dz 


^ + k^^(z)  - k^^  COS^  0 - -J 
- dz^^ 


= 0 , 


with  k^(z)=  cn/c^(z)  and  From  Eq.  (!>)  it  follows  that 


(z) 


^ . ^1' 

2 x^^  (z) 


(5) 


(6) 


Using  Eq.  (6)  in  Eq.  (4)  the  reflection  coefficient  may  be  expressed  as 

R = ” ^ |z=u  (7) 
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with  a = p /p^(0)  and  K = k sin  6. 
o 1 o o 

A first  order  treatment  of  Eq.  (7)  consists  of  neglecting  the 
effect  of  the  density  variation  on  the  velocity  profile,  that  is, 


av 

IK  - iiL  _ a — 

o 2 


R = 


(0) 


z=0 


i/C  + ^ + a 


X;;^’(0) 

Xi(o)'' 


z=0 


where  x^(0)(z)  satisfies 


^ + k ^(z)  - k ^ cos^  B 
,2  1 o 

dz 


\(0)(z)  = 0 


The  expression  for  R may  now  be  written  as 


N 

R =-°  " 


D 

o z 


where 


N = iK  - a X,  ' ( 0 ) / Xt  ( 0 ) I 
o o 1 1 


z=0 


The  expression  for  R may  now  be  written  as 
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N “ 2N  •' 

R = = R 
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gr 

2N 
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'o[l.|f] 
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where  R is  the  reflection  coefficient  when  p ' (z)=0. 
o 1 

Since  T is  in  general  very  small,  it  is  possible  to  approximate  the  squared 
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magnitude  of  R as  follows 


Equation  (lO)  indicates  that  the  direction  of  the  shift  in  jR|  due  to 
a density  gradient  is  controlled  by  the  sign  of  the  phase  of  the  reflection 
coefficient  for  p'(z)=0.  To  see  this  effect,  consider  a 100  ra  clay  layer 
with  a constant  sound  speed  gradient  of  1 sec  ^ overlying  a rock  sub- 
strate. The  parameters  for  the  rock  and  clay  are  listed  in  Table  II-2. 
Figure  11-26  depicts  bottom  loss  curves  for  this  clay  layer  for  three 
different  values  of  constant  densitr  gradient,  p'=(0,  0.001,  and  0.002) 
(g/cm^)/m.  Since  only  rays  encountering  the  bottom  at  grazing  angles 
greater  than  18.78°  will  encounter  the  substrate,  the  above  derived  results 
should  be  valid  in  the  small  grazing  angle  region.  Figure  11-27  depicts 
the  typical  phase  behavior  of  the  reflection  coefficient  for  a 100  m 
clay  layer  for  p'=0.  This  figure  shows  that  the  phase  approaches  -n 
through  the  third  quadrant  as  d goes  to  zero.  Hence,  Eq.  (lO)  indicates 
that 


jR|^  < |Rq|^  as  e 0 

Since  bottom  loss  is  defined  as  BL  = -20  log  |r|,  it  follows  that 


bl(|r| ) > BL( |R^| ) as  e -»  0 

This  is  precisely  the  behavior  of  Fig.  11-26  for  small  6.  Figures  11-28 
and  11-29  depict  bottom  loss  curves  for  silt  and  sand  layers  overlying 
rock  substrate.  The  phase  of  R for  these  sediments  behave  similarly  to 
that  of  clayj  hence  as  9 0,  BL(r)  > BL(R^).  Although  Figs.  11-28  and 

11-29  do  not  clearly  show  it,  this  is  actually  the  case.  The  differences 
are  so  small,  however,  that  they  could  not  be  resolved  on  the  plotter. 
Needless  to  say,  in  the  absence  of  subbottom  reflectors,  the  effects  of  a 
density  gradient  are  negligibly  small. 

3 . Boundary  Condition  Effects 

It  is  important  to  remember  that  the  effects  just  described  are 
due  to  the  changes  in  form  of  the  wave  equation  itself,  i.e.,  the  presence 
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BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  A 100  tn  CLAY 
SEDIMENT  LAYER  OVERLYING  A ROCK  SUBSTRATE 


FREQUENCY:  50  H 


FIGURE  ir-27 

REFLECTION  COEFFICIENT  PHASE  versus  GRAZING  ANGLE 
FOR  A SOLID  SUBSTRATE 
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FIGURE  n-29 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  A 100  m SAND 
SEDIMENT  LAYER  OVERLYING  A ROCK  SUBSTRATE 


of  terms  in  the  wave  equation  involving  the  density  gradient.  The  indirect 
effect  due  to  altered  boundary  conditions  is  to  be  considered  next. 

How  a density  gradient  affects  the  bottom  loss  at  larger 
grazing  angles  may  be  understood  by  considering  the  changes  of  Impedance 
at  subbottom  reflecting  horizons,  specifically  the  rock  substrate  in  the 
cases  considered  here.  In  general,  since  rock  density  is  usually  higher 
than  sediment  densities,  the  effect  of  a positive  density  gradient  is 
to  increase  the  impedance  matching  across  a sediment-rock  interface.  This 
increased  impedance  matching,  in  turn,  enhances  the  transmission  into  the 
substrate  which,  in  turn,  increases  the  loss.  This  effect  will  begin 
when  waves  impinge  upon  the  substrate  at  near  the  shear  wave  critical 
angle  since,  in  this  region,  transmission  into  the  substrate  becomes  non- 
negliglble. 

For  the  clay,  silt,  and  sand  examples  being  considered,  the 
bottom  angles  corresponding  to  the  shear  wave  critical  angle  and  the 
compressional  wave  critical  angle  are  55-22°  and  7^.55° ^ respectively. 
Figures  11-26  and  11-28  show  the  expected  results,  i.e.,  higher  bottom 
loss  in  the  presence  of  density  gradients  above  55-22°. 

Figure  11-29  for  sand  shows  the  same  general  features  above  the 
shear  wave  critical  angle  with  the  exception  of  the  peaked  regions  at 
about  60°  and  78°.  In  these  regions  the  bottom  loss  for  p'ijio  is  less  than 
for  p'=0.  Sand  is  different  from  clay  and  silt  in  that  it  has  higher 
attenuation  and  it  is  a high  speed  bottom  exhibiting  a critical  angle 
effect.  It  seems  possible  that  some  coupling  between  the  attenuation 
and  density  gradient  might  be  responsible  for  the  deviation  around  78°. 

Figure  11-50  depicts  the  bottom  loss  for  sand  using  an  attenuation 
reduced  by  an  order  of  magnitude.  With  the  exception  of  the  region  around 
60°,  the  effect  has  been  reversed  and  is  now  similar  to  a silt  or  clay. 

The  interplay  between  p'  and  attenuation  and  the  nature  of  the  reversal 
effect  at  60°  are  topics  for  further  consideration. 
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FIGURE  11-30 

BOTTOM  LOSS  versus  CRAZING  ANGLE  FOR  A 100  m,  LOW  ATTENUATION 
SAND  SEDIMENT  LAYER  OVERLYING  A ROCK  SUBSTRATE 


To  test  the  hypothesis  that  the  density  gradient  can  be  handled 
to  a good  approximation  by  considering  only  boundary  condition  effects, 
the  clay  bottom  loss  was  redetermined  neglecting  the  effect  of  the  density 
gradient  on  the  wave  equation.  In  other  words,  the  pressure  variations 
throughout  the  clay  layer  were  calculated  niamerically  as  if  p'=0  whereas 
the  boundary  conditions  to  be  satisfied  were  specified  assuming  a nonzero 
density  gradient.  The  bottom  loss  computed  with  this  pressure  function 
was  then  compared  to  that  computed  with  the  pressure  field  calculated 
with  a full  treatment  of  the  density  gradient.  The  results  of  this 
comparison  for  p' =0.002  (g/cm^)/m  were  excellent.  In  the  low  angle  region 
(0°  to  10°)  where  boundary  condition  effects  are  expected  to  be  unimportant, 
the  most  deviation  was  noted;  the  average  deviation  in  the  bottom  loss  was 
0.05  ±0.05  dB.  In  the  higher  angle  regions  (20°  to  90°)  where  the  boundary 
effects  are  important,  the  differences  were  less  and  were  too  small  to  be 
graphically  displayed.  The  average  deviation  in  bottom  loss  was 
0.01  ± 0.008  dB.  The  angular  region  from  10°  to  20°  was  excluded  because 
of  the  anomalous  behavior  in  this  region.  (The  large  peak  at  about  1^° 
is  due  to  the  excitation  of  a Stoneley  wave.) 

Several  dominant  features  are  apparent  in  the  previously 
mentioned  bottom  loss  curves,  e.g.,  the  prominent  peaks  at  low  angles 
for  silt  and  clay  and  the  peak  at  about  53°  in  clay.  The  peak  in  clay 
bottom  loss  at  55°  is  associated  with  the  excitation  of  a Rayleigh  wave. 

A Rayleigh  wave  is  an  Interface  wave  which  exists  Ideally  at  the  free 
surface  of  a solid.  The  enhancement  of  this  peak  when  p'=0  is  consistent, 
since,  in  this  situation,  the  sediment  density  above  the  substrate  is 
lower  and  the  substrate  boundary  more  closely  approaches  a free  surface. 

This  Rayleigh  wave  peak  is  also  present  in  silt  although  its  magnitude 
does  not  show  the  sensitivity  to  p'  that  clay  does.  This  is  likely  related 
to  the  fact  that  silt  attenuation  is  much  higher  than  it  is  in  clay. 

The  low  angle  peak  at  about  l4°  in  clay  has  been  previously 
identified  as  due  to  excitation  of  a Stoneley  wave  (Hawker^) . A Stoneley 
wave  is  an  Interface  wave  excited  on  the  surface  of  a solid.  Figure  11-26 
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shows  that  a density  gradient  tends  to  shift  the  location  of  the  Stoneley 
wave  peak.  The  ahovementioned  procedure  of  using  only  interface  densities 
to  treat  a density  gradient  was  able  to  account  for  the  shift  in  peak 
location.  It  is  expected  that  further  work  relating  this  shift  to  density 
gradients  and  other  sediment  parameters  will  be  done. 

The  bottom  loss  curve  for  silt  shown  in  Fig.  11-28  has  a feature 
which  resembles  a Stoneley  wave  peak  at  about  l4°.  This  feature  is  not, 
however,  due  to  excitation  of  a Stoneley  wave  since  the  feature  remains 
intact  when  the  substrate  is  treated  as  a liquid.  This  feature  has  an 
interesting  behavior  in  that  its  prominence  decreases  with  increasing  p', 
whereas  its  location  remains  unchanged.  The  nature  of  this  feature  in 
silt  is  a topic  for  further  consideration  in  future  research. 

4.  Conclusions 


A STirvey  of  Figs.  11-26,  11-28,  and  11-29  indicates  that,  in 
general,  the  density  gradient  affects  clay  the  most,  silt  to  a lesser 
degree,  and  sand  to  the  least  degree.  With  the  exception  of  the  low  angle 
(=14°)  anomalous  behavior  of  silt  and  clay,  the  nonnegligible  effects  are 
in  the  higher  grazing  angle  regions.  Since  the  effects  in  this  region  are 
"boundary  effects,"  it  seems  likely  that  the  relative  magnitude  of  the 
density  gradient  effect  should  be  related  to  the  fractional  change  of 
impedance  at  the  substrate  interface  when  density  gradients  are  present 
relative  to  the  impedance  in  the  absence  of  a density  gradient. 

Table  II-3  shows  the  values  of  the  sediment  characteristic 
Impedances  just  above  the  rock  substrate  for  100  m clay,  silt,  and  sand 
layers  having  sound  velocity  gradients  of  1 sec  The  numbers  in  paren- 
theses are  the  percent  differences  of  the  impedances  for  p' =0.0001  and 
0.002  relative  to  p'=0.  The  trend  of  these  ratios  indicates  a larger 
relative  effect  for  clay,  followed  by  silt,  and  then  sand.  This  trend  is 
borne  out  in  the  bottom  loss  model  results.  For  p'=0.002  (g/cm  ,/“>  4he 
percent  differences  in  bottom  loss  for  the  clay,  silt,  and  sand  layers  in 
the  angular  range  of  55°  4o  90°  relative  to  p'=0  were  12.5^>  7-0^,  and 
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TAB1£  II-5 


SEDIMENT  CHARACTERISTICS  IMPEDANCES  EVALUATED  JUST 
ABOVE  THE  SUBSTRATE 

2 

Sand  kg/(m  -sec) 

p'  = 0 (g/cm^)/ni  '>-79  X 10^ 

p'  = 0.001  (g/cm^)/m  5-98  X 10^ 

p'  = 0.002  (g/cni^)/ra  4.l6  X 10^  (9-6^) 

sediment  sound  speed  gradient  = 1 sec”^ 

sediment  layer  thickness  = 100  m 


2 2 > 
Slit  kg/(m  -sec)  Clay  kg/(m  -sec! 


.72  X 10 


, 6 
2.06  X 10 


5.88  X 10^  (5.9^)  2.25  X 10^  (7-5^) 


5.0A  X 10^  (11. 9'^)  2.59  X 10^  (15.7^) 
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respectively.  While  the  absolute  validity  of  this  procedure  might  be  suspect, 
especially  in  the  light  of  a possible  coupling  between  p'  and  other  subbottom 
parameters,  it  does  seem  to  provide  a simple  method  of  gauging  the  density 
gradient  effects. 

In  the  low  angle  region,  which  is  actually  of  m'^re  Interest  in 
long  range  propagation  studies,  there  is  little  that  can  be  said  about 
density  gradient  effects  on  a sediment-by-sediment  basis.  In  this  region 
the  effects  are  usually  buried  in  the  wave  equation  and  are  not  easily 
determined  except  near  9=0° . 

It  is  worth  emphasizing  at  this  point  that  the  density  gradient 
effects,  whether  in  the  high  or  low  angle  regions,  are  small  effects.  In 
fact,  one  goal  of  this  and  future  research  is  to  determine  whether  it 
is  ever  necessary  to  account  for  the  density  gradient,  i.e.,  is  the 
effect  of  a density  gradient  ever  more  than  negligible?  Since  only  three 
sediment  types  and  one  layering  configuration  were  studied  here,  it  is  too 
soon  to  draw  firm  conclusions.  However,  based  on  the  results  quoted  here, 
it  is  possible  to  make  a few  observations. 

In  the  high  angle  region,  the  variations  due  to  introduction  of 
a density  gradient,  in  general,  occur  over  a wider  angular  range,  as 
opposed  to  variations  in  narrow  peaked  regions  in  the  low  angle  range. 

Hence,  the  potential  for  causing  noticeable  effects  on  propagation  in  the 
water  column  are  the  greatest.  It  seems  likely,  however,  that  there  might 
be  other  effects  which  tend  to  counter  this  potential;  one  such  effect 
is  sediment  thickness.  As  the  sediment  layer  becomes  thicker,  it  is 
conceivable  that  the  added  attenuation  loss  might  tend  to  reduce  the 
relative  density  gradient  effect. 

Another  factor  which  tends  to  lessen  the  effect  of  density 
gradients  in  the  high  angle  region  is  the  diminishing  Importance  of  the 
high  grazing  angles  in  long  range  propagation  studies. 

In  the  low  angle  region  (less  than  20°)  the  examples  considered 
in  this  study  tend  to  show  that  the  only  time  a density  gradient  might  be 
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significant  is  when  it  is  associated  with  some  sort  of  anomalous  behavior 
such  as  interface  waves,  intromission  angles,  etc.  Since  these  anomalies 
tend  to  be  fairly  narrow  in  width,  it  is  uncertain  as  to  whether  the 
effects  due  to  density  gradients  would  ever  be  noticeable  at  these  angles. 

D.  Additional  Bottom  Loss  Studies 

In  addition  to  the  studies  reported  in  the  first  three  sections  of 

this  chapter,  various  associated  studies  were  carried  out  during  this 

contract  year.  This  section  summarizes  briefly  some  of  the  principal 

results  found  in  the  studies  of  hidden  depth  and  the  sensitivity  of  bottom 

loss  variations  in  sediment  type.  Both  of  these  subjects  have  been  report- 

2 

ed  in  detail  by  Hawker,  Focke,  and  Anderson. 

1.  Hidden  Depth  Studies 

Although  the  concept  of  a hidden  depth  intuitively  appears  to 
be  clear,  a more  precise  and  quantitative  definition  will  be  required  if 
specific  results  are  to  be  obtained. 

One  approach  to  the  problem  is  via  mode  theory,  with  the  hidden 

depth  defined  as  the  depth  below  which  the  fractional  change  in  the 

eigenvalues,  caused  by  the  introduction  of  a perfect  reflector,  is  less 

than  a specified  tolerance.  This  approach  necessitates  specifying  which 

modes  are  to  be  considered  and  therefore  is  not  independent  of  the  pro- 

12 

pagatlon  geometry.  This  approach  is  the  one  taken  by  Williams  In  his 
investigation  of  hidden  depths. 

An  alternative  approach,  and  the  one  adopted  in  this  study,  is 
to  define  the  hidden  depth  in  terms  of  the  fractional  change  in  the  bottom 
loss  caused  by  the  introduction  of  a perfect  reflector  at  some  depth. 

When  this  change  is  less  than  a certain  tolerance,  the  hidden  depth  has 
been  reached.  Here  again  some  elements  of  the  propagation  geometry  are 
hidden  since  the  grazing  angle,  or  angular  interval,  must  be  specified. 
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Also,  in  addition  to  the  freedom  to  choose  the  tolerance,  there 
is  the  possibility  of  modifying  the  definition  to  examine  the  difference 
between  a normal  sediment  colvimn  and  one  with  an  absorber  underneath,  or 
perhaps  the  difference  between  the  reflector  and  absorber  cases,  'fhe 
most  extreme  contrast  is  obtained  by  comparing  situations  in  which  a 
column  of  sediment  overlays  alternately  a perfect  reflector  and  a perfect 
absorber.  This  definition  will  be  adopted  in  this  report.  That  is,  the 
hidden  depth,  Zj^,  is  defined  as  the  depth  at  which  a reflector  and  then 
an  absorber  can  be  inserted  in  the  sediment  such  that  = 1 20  log 

/|r|a-u  a-  where  IrI  is  the  modulus  of  the  plane 

''  'Reflection  ' 'Absorber'^'  ' ' 

wave  reflection  coefficient  and  € is  the  tolerance.  Investigation  has 
shown  that  the  precise  nature  of  the  comparison  cases  is  not  of  great 
importance  and  the  results  would  be  little  changed  if  the  definition  were 
somewhat  altered. 

In  the  remainder  of  this  subsection  we  shall  examine  the  hidden 
depth  problem  for  two  sediment  types,  a fine  clay  and  a coarse  silt.  The 
clay  sediment  has  a density  ratio  (relative  to  sea  water)  of  1.206,  a sound 
speed  ratio  of  0.991^  and  a porosity  of  87^.  The  silt  sediment  has  a 
density  ratio  of  1.598,  a sound  speed  ratio  of  O.986,  and  a porosity  of 
65^.  The  frequency  was  chosen  to  be  50  Hz;  at  this  frequency  the  attenua- 
tions were  O.OO285  ds/m  for  clay  and  0.012  ds/m  for  silt.  The  sediment 

properties  and  their  relationship  to  porosity  follows  the  work  of 
7 

Hamilton. 


Flgiires  II-5I  through  11-55  show  the  bottom  loss  versus  grazing 
angle  versus  layer  thickness  for  a clay  layer  overlying  a perfectly  re- 
flecting surface.  These  bottom  loss  curries,  as  well  as  all  others  appear- 
ing in  this  report,  were  computed  using  the  bottom  loss  model  developed  at 
ARLrlTT  by  Hawker  and  Foreman.^  The  sound  speed  profile  was  assumed  to  be 
linear  C(z)=C(0)+gz,  and  the  gradients  g are  0.5  sec  1.0  sec  and 

1.5  sec”^  in  the  three  figures.  This  range  of  sound  speed  gradients  in- 
cludes those  most  commonly  found  in  deep  ocean  sediments. 
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BOTTOM  LOSS  - dB 


BOTTOM  LOSS  - dB 


These  three  figures  all  display  a prominent  diagonal  ridge  with 
respect  to  the  angle  axis'  decreasing  with  increasing  gradient.  For  any 
given  angle,  at  depths  beyond  this  ridge,  the  bottom  loss  is  seen  to  be 
independent  of  any  further  increase  in  layer  thickness.  This  indepen- 
dence of  layer  thickness  is  an  indication  that,  at  this  angle,  the  hidden 
depth  has  been  reached.  Moreover,  the  decreasing  inclination  of  the  ridge 
shows  that,  at  a fixed  angle,  the  hidden  depth  is  reached  more  quickly  for 
greater  sound  speed  gradients.  This  effect  is  simply  due  to  the  increased 
upward  refraction  of  the  soiind  speed  profile  for  larger  gradients. 

Consider  now  a quantity,  A(6,z),  defined  as  the  absolute  value 
of  the  difference  between  the  bottom  losses  in  the  case  of  a perfectly 
reflecting  plane  and  a perfect  absorber  inserted  at  a depth  z.  This 
difference  can  be  seen  to  be  simply 

= I -20  log  (I  I Reflector^  ^ Absorber^  ^ 

= 20 1 log  (iRlpefigctor/I^Ubsorber^  I 

which  is  simply  the  intensity  ratio  in  the  two  cases  expressed  in  decibels. 

The  hidden  depth,  defined  by  A(0,z)^e,  is  directly  obtainable  from 
A(0,z)  once  the  tolerance  e is  chosen.  Since  there  is  no  particular  rea- 
son to  choose  € to  have  any  particular  value  in  general  (in  a specific 
problem  it  might  be  fixed),  it  might  seem  at  first  that  a considerable 
element  of  arbitrariness  remains  in  the  problem.  However,  it  has  been 
shown  that  the  s\irface  A(e,z)  is  displayed  as  a sharp  transition  in  the 
vicinity  of  the  hidden  depth.  As  a result  of  this  steep  transition,  a 
change  of  e from,  say,  0.1  dB  to  0.01  dB,  will  cause  very  little  change  in 
the  derived  value  of  the  hidden  depth.  In  other  words,  the  concept  of  a 
hidden  depth  will  be  a useful  one  and  the  tolerance  e need  not  be  specified 
very  precisely. 

Figure  shows  the  hidden  depth,  defined  with  e = 0.1  dB, 
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FIGURE  ir-34 

HIDDEN  DEPTH  versus  GRAZING  ANGLE  FOR  A CLAY 
LAYER  WITH  VARIOUS  SOUND  SPEED  GRADIENTS 


plotted  versus  grazing  angle  for  the  three  sound  speed  gradients  shovm  in 
the  preceding  groups  of  figures.  For  comparison  purposes.  Fig.  11-55 
shows  the  ray  txirning  depth  versus  grazing  angle  for  these  same  gradients. 
It  will  be  observed  that  the  hidden  depth  always  lies  several  wavelengths 
below  the  turning  depth.  This  close  association  of  hidden  depth  and  turn- 
ing depth  is  not  surprising  since  it  is  at  the  turning  depth  that  the 

pressure  field  changes  from  oscillatory  to  exponentially  decaying  behavior. 

12 

In  the  investigation  carried  out  by  Williams,  the  hidden  depth  was 
in  fact  referenced  to  the  turning  depth,  and  was  again  found  to  lie  a 
few  wavelengths  below  it. 

Having  established  the  hidden  depth  in  one  particular  case, 
we  come  to  the  general  question  of  the  variation  of  the  hidden  depth 
with  changes  in  sediment  type,  with  different  sound  speed  profiles, 
with  frequency,  etc.  Althou^  these  questions  have  not  yet  been  completely 
resolved  and  these  lines  of  investigation  are  still  being  pursued,  some 

t 

results  have  been  obtained  and  are  reported  by  Hawker,  Focke,  and  Anderson. 
A single  such  example  is  contained  in  Fig.  II-36  which  shows  the  hidden 
depth  for  a silt  layer  for  three  sound  speed  gradients.  These  results 
are  generally  similar  to  those  for  clay  shown  in  Fig.  11-3^. 

The  close  association  between  the  hidden  and  turning  depths 
leads  at  once  to  the  conclusion  that  the  apparently  very  complex  dependence 
of  hidden  depth  on  all  of  the  various  sediment  parameters  is  in  actuality 
not  an  overriding  problem.  Upon  writing  the  hidden  depth  as  'Z^= 
where  k is  some  sound  wavelength  (say  the  water  wavelength)  and  6 Is  a 
parameter,  we  see  that  all  of  the  complexity  is  contained  in  & alone. 
Moreover,  we  now  have  reason  to  believe  that  6 is  of  order  unity,  that  is, 
it  is  not  much  larger  than  one.  In  general  & will  depend  upon  all  the 
sediment  parameters  as  well  as  frequency.  The  dominant  sound  speed  pro- 
file dependence  should  be  contained  in  the  term  Z^. 

Thus  far  we  have  restricted  our  investigation  to  the  case  of  a 
linear  sound  speed  profile.  Althou^  it  is  true  that,  in  the  uppermost 
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FIGURE  ir-35 

TURNING  DEPTH  versus  GRAZING  ANGLE  FOR  A CLAY 
LAYER  HAVING  VARIOUS  SOUND  SPEED  GRADIENTS 
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sediment  layers,  the  so\md  speed  gradient  is  generally  constant,  seismic 
profiling  measurements  clearly  show  that  the  gradient  decreases  with  depth. 

To  gain  some  insight  into  how  such  behavior  might  affect  the  hidden  depth 
picture,  we  shall  examine  the  exponential  sound  speed  profile  proposed 
by  V/illiams,  given  by  C(z)=C(0)[l+A  (e  ^^-l)]  This  profile  is  a 

3-parameter  one  with  C(<»)=C(0) [ 1-A^]  C'(Z)=C'(o)  e ^^[l+A^(e  ^^-l)]  , 

and  C (o)=l/23A^C(o) . Figure  11-37  shows  C'(Z)  plotted  versus  depth  for 
various  values  of  C(o)/c(a)),  with  C'(0)  chosen  to  be  unity.  It  will  be 
observed  that  for  A^<i2/3',  C'(z)<l  for  all  depths.  Figure  II-38  shows  the 
ray  turning  depth  for  the  case  of  water  overlying  sediment  containing 
such  an  exponential  sound  speed  profile.  As  would  be  expected  from  Fig. 

II-38  for  C(0)/c(<o)>l/»^,  the  ray  turning  depth  will  always  be  larger  than 
it  is  in  the  case  of  a linear  profile  (dotted  line).  For  a given  value 
of  C(0)/C(a>),  the  turning  depth  approaches  infinity  at  a certain  angle. 

These  critical  angles  are  denoted  by  the  asymptotes. 

The  most  important  result  contained  in  these  figures  concerns 
the  deviation  of  Z^  from  linear  behavior  at  low  angles.  It  will  be  noted 
that,  for  e§25°  and  0.85C(0)/C(«>)&  1/*^  the  maximum  deviation  from  linear 
behavior  is  ^5  m.  For  this  range  of  C(0)/c(“),  the  sound  speed  gradient 
at  200  m can  range  from  0.66  to  0.9T!><  Thus,  although  a considerable 
variation  of  C'(z)  can  occur  the  turning  depth  will  range  only  over  a 
comparatively  narrow  Interval.  Consequently,  the  hidden  depth  values 
which  are  previously  determined  in  several  cases  using  the  linear  profile 
will  not  be  greatly  altered  by  the  introduction  of  such  an  exponential 
profile . 


2.  Sensitivity  of  Bottom  Loss  to  Variations  in  Sediment  Type 

Even  with  all  other  problem  parameters  fixed  (layering,  sound 
speed  profile  function  and  gradient,  substrate,  etc.)  the  problem  of  the 
influence  of  sediment  type  variations  on  bottom  loss  Involves,  at  a mini- 
mum, variations  in  sound  speed,  density,  and  attenuation.  One  possibility 
for  circumventing  this  difficulty  is  to  relate  these  three  parameters  to 
a single  parameter  such  as  porosity  or  mean  grain  size.  The  work  of 
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Hamilton  and  Akal  has  shown  that,  for  a considerable  variety  of  sediments, 
these  three  parameters  can  be  usefully  related  to  porosity  and  the  appro- 
priate equations  have  been  developed  by  these  authors. 

Upon  fixing  all  other  parameters  in  the  problem  then,  it  is 
possible  to  study  sediment  type  variations  be  examining  the  behavior  of 
bottom  less  versus  angle  versus  porosity.  An  empirical  approach  to  this 
type  of  study  was  taken  by  Hall  and  Watson^^  who  used  empirical  fits  to 
AMDS  data  and  developed  curves  giving  bottom  loss  versus  grazing  angle 
parameterized  by  porosity. 

It  is  quite  straightforward  to  employ  the  aforementioned  results 
of  Hamilton  and  Akal  to  generate  appropriate  parameter  sets  for  use  in  our 
bottom  loss  model.  This  has  been  done  and  Figs.  11-59  and  II-4o  show  two 
resulting  3-<ilmensional  surfaces  in  the  case  of  a linear  sound  speed  pro- 
file with  a gradient  of  1.2  sec  ^ in  a 100  and  200  m layer,  respectively. 

The  substrate  is  a homogeneous  solid. 

There  are  two  striking  features  of  these  surfaces,  the  bladed 
structures  at  low  angles  and  the  c\irvlng  ridge  at  higher  angles.  The 
prominent  curving  ridge  is  associated  with  Stoneley  waves  at  the  sediment- 
substrate  interface  whereas  the  bladed  structures  seem  to  be  a residue 
of  an  Intromission  angle  effect.  It  is  noteworthy  that,  since  the  bladed 
structures  lie  at  angles  such  that  the  substrate  is  below  the  hidden 
depth,  no  change  is  seen  in  these  structures  between  the  100  m and  200  m 
cases . 


These  curves  show  directly  how  the  bottom  loss  can  be  expected 
to  change  with  sediment  type  variations.  Future  work  will  include  pro- 
ducing additional  such  curves  with  different  sound  speed  gradients  and 
at  different  frequencies. 
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FREQUENCY:  50  Hi 

LAYER  THICKNESS:  200  m 

SOUND  SPEED  GRADIENT:  1.2  sec"’ 


FIGURE  ir-40 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
versus  POROSITY  FOR  A 200  m LAYER 
OVERLYING  A SOLID  SUBSTRATE 
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III.  SENSITIVITY  OF  SOUND  PROPAGATION  TO  VARIATION 
IN  BOTTOM  LOSS 


The  objectives  of  this  study  allow  the  bottom  interaction  problem 
to  be  defined  by  the  questions  "what  influence  does  the  bottom  have  on 
long  range  acoustic  propagation?"  and  "how  much  information  about  the 
(sub)bottom  is  required  in  order  to  adequately  predict  these  effects?" 
This  restricted  form  of  the  problem  is  still  quite  complex  and  requires 
the  problem  to  be  separated  into  simpler  components. 

A systematic  sensitivity  study,  summarized  in  this  chapter,  was 
conducted  to  isolate  regions  in  the  water  column  exhibiting  strong  bottom 
influence  from  regions  which  are  relatively  insensitive  to  bottom  effects. 
For  this  purpose  simple  bottom  loss  descriptions  were  used  to  control  the 
bottom  interactions.  Throughout  most  of  this  study  variations  in  the 
bottom  loss  were  restricted  to  changes  in  critical  angle:  bottom  loss 
was  0 dB  below  the  critical  angle  and  essentially  infinite  above.  Many 
of  the  results  here  have  been  given  by  Hawker,  Focke,  and  Anderson.^ 

A . Methodology 


Two  propagation  models  were  used  in  this  investigation,  FACT,  a ray 

2 3 

trace  model  reported  by  Spofford  and  Baker  and  Spofford,  and  a parabolic 

4 

equation  model  reported  by  Brock.  The  nature  of  the  ARL:UT  implemen- 

5 

tation  of  these  models  is  described  by  Hawker,  Foreman,  and  Focke.  FACT 
defines  the  bottom  in  terms  of  a plane  wave  reflection  coefficient  applied 
at  the  water- sediment  Interface.  The  parabolic  equation  model  simulates 
a critical  angle  effect  by  introducing  below  the  water  column  a thin  layer 
with  a sound  speed  gradient  and  zero  attenuation,  which  in  turn  overlies 
a homogeneous  attenuating  layer. 


Various  terras  used  throughout  this  chapter  are  defined  here.  The 


depth  at  which  the  sound  speed  is  the  same 
referred  to  as  the  source  conjugate  depth. 
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as  that  at  the  soxirce  depth  is 
The  depth  at  which  the  sound 
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speed  is  the  same  as  that  at  the  surface  is  called  the  critical  depth. 

The  distance  between  the  critical  depth  and  the  bottom  (if  the  bottom  lies 
below  critical  depth)  is  called  the  depth  excess. 

A sensitivity  study  is  designed  to  measure  the  influence  of  a given 
input  parameter  on  a specified  output  parameter.  For  the  study  discussed 
in  this  chapter,  systematic  variations  are  introduced  in  the  bottom  loss 
to  determine  their  effects  upon  the  total  calculated  acoustic  propagation. 
The  major  Intensity  contribution  from  bottom  interacting  energy  will  be 
within  the  shadow  zones,  away  from  the  convergence  zones  where  waterborne 
paths  dominate  the  intensity  calculations.  Figure  III-l  demonstrates 
this  bottom  influence  on  the  energy  contributions  within  the  shadow  zones. 
This  range  dependence  of  the  influence  of  the  bottom  (dominant  only  in  the 
shadow  zones)  is  a function  of  the  environment  and  the  source/receiver 
geometry.  To  remove  problems  associated  with  this  range  dependence  in 
a manner  independent  of  the  environment,  average  propagation  losses  were 
computed  as  the  output  of  primary  importance.  Major  conclusions  from  the 
sensitivity  study  are  based  on  the  resultant  variations  in  these  averages. 

The  form  of  averaging  used  in  this  study  computes  the  propagation  loss 
for  the  average  intensity  at  a given  receiver  depth  within  the  100  to  200 
nmi  range  interval.  This  range  Interval  will  not  include  the  nearfield 
effects  (the  direct  ray  paths).  It  will  cover  several  convergence  zones 
and  will  have  only  a nominal  3 dB  increase  in  loss  due  to  increases  in 
range.  Run  times  for  the  P.E.  model  will  not  be  excessive  out  to  200 
nmi. 


Variations  in  the  bottom  loss  were  limited  to  changes  in  the  critical 
angle,  using  values  of  0°,  1°,  5°/  10°/  15°/  and  20°.  Figure  III-l  pre- 
sents the  propagation  loss  versus  range  using  these  various  bottom  descrip- 
tions for  a given  soijrce/receiver  depth  combination.  The  profile  used  is 
from  the  mid- Pacific  Ocean.  It  has  a critical  depth  of  3952  m,  an  axis 
depth  of  600  m,  and  a 44  m/sec  difference  in  sound  speed  between  the  surface 
and  the  axis.  As  expected,  the  bottom  descriptions  with  the  higher  critical 
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FIGURE  III-l 

PROPAGATION  LOSS  FOR  VARIOUS  BOHOM  LOSS  CRITICAL  ANGLES 
SOURCE  DEPTH:  152  m,  RECEIVER  DEPTH:  3350  m,  BOHOM  DEPTH:  3952  m.  MODEL:  FACT 


angles  contribute  more  energy  within  the  shadow  zones.  The  maximum 
acoustic  intensities  within  the  convergence  zones,  however,  remain  virtually 
unchanged  for  the  bottom  variations. 

B.  Results 

Average  propagation  losses  were  computed  at  ninnerous  receiver  depths 
for  a selected  environment  and  a given  source  depth.  These  averages  were 
then  plotted  as  a function  of  receiver  depth.  These  average  propagation 
loss  cxirves  computed  by  FACT  are  presented  in  Fig.  III-2  for  the  same 
environment  as  in  Fig.  III-l.  Again,  as  expected,  an  increase  in  critical 
angle  decreases  the  average  propagation  loss. 

A source  depth  of  152  m has  a conjugate  depth  of  2500  m.  Each  of  the 
curves  in  Fig.  III-2  has  a relative  minimum  propagation  loss  at  both  these 
depths  (resulting  from  cusped  caustice.  Ref.  2).  Outside  the  region  be- 
tween these  two  depths,  the  average  propagation  loss  increases  as  the  re- 
ceiver depth  approaches  either  the  surface  or  the  bottom. 

The  intensity  contributions  from  the  various  bottom  types  can  be 
more  readily  observed  by  taking  the  propagation  loss  differences  between 
the  0°  critical  angle  case  and  each  of  the  remaining  bottom  type  cases. 
Figure  III-3  presents  these  differences  for  the  data  previously  seen  in 
Fig.  III-2.  From  Fig.  III-3  the  water  column  appears  to  be  divided  into 
three  separate  regions:  from  the  surface  to  the  source  depth  (l52  m), 
between  the  source  depth  and  Its  conjugate  depth  (2500  m),  and  from  the 
sotirce  conjugate  to  the  bottom  (3952  m) . Within  the  first  and  third 
regions  the  bottom  Influence  Increases  as  the  receiver  depth  approaches 
either  the  stirface  or  the  bottom.  The  bottom  loss  contributions  increase 
as  the  receiver  depth  approaches  either  boundary  for  each  of  the  specified 
bottom  types  (an  anomaly  does  occur  at  the  bottom  for  the  0“  critical 
angle  case,  creating  maximum  differences  at  the  bottom  depth). 

Within  the  second  region  of  the  water  coliunn,  the-  bottom  influence 
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FIGURE  III-2 
AVERAGE  PROPAGATION  LOSS 

Range  Interval:  100-200  nm 


is  almost  independent  of  the  receiver  depth.  The  propagation  loss  has  a 
regular  dependence  on  the  bottom  loss  within  this  region.  Table  III-l 
presents  a comparison  between  the  propagation  differences  in  this  second 
region  and  the  bottom  bounce  energy.  This  bottom  bounce  energy  is  the 
additional  energy  Introduced  into  the  problem  by  the  bottom  bounce  rays 
in  the  form  of  total  energy  divided  by  water  refracted  energy.  These 
niombers  are  based  on  Snell's  law  using  the  sound  speed  at  the  source 
depth  and  at  the  bottom.  Within  this  second  region  of  the  water  coliamn, 
there  is  at  most  a 0.2  dB  difference  between  the  intensity  contribution 
and  the  energy  contributions  for  the  bottom  bounce  rays  in  this  one  case. 

Similar  calculations  were  performed  for  various  source  depths. 

Figure  III-4  presents  the  proj«gation  loss  differences  between  the  0°  and 
20“  critical  angle  cases  for  five  source  depths.  The  conjugate  depths  for 
these  five  source  depths  are  387^  m (55*5  3^10  (50  2963  m (90  m) , 

2500  m (152  m),  and  600  m (6OO  m).  For  each  source  depth  the  water  column 
is  divided  into  the  regions  defined  by  the  soiirce  depth  and  its  conjugate 
depth.  The  characteristics  described  earlier  are  seen  to  be  consistent 
for  each  of  these  source  depths. 

Table  III-2  presents  comparisons  of  the  bottom  bounce  energy  and  the 
bottom  bounce  contributions  to  the  average  propagation  loss.  The  two 
calculations  differ  by  less  than  1.0  dB,  except  at  the  axis  depth  of 
600  m.  (The  6OO  m source  is  a sp)eclal  case  where  the  source  depth  is  its 
own  conjugate  depth.)  The  bottom  bounce  contributions  to  the  average 
propagation  loss  can  be  estimated  by  a calculation  of  the  bottom  bounce 
energy  for  receiver  depths  between  the  source  and  source  conjugate  depths. 

Additional  studies  have  been  conducted  using  this  same  mld-Pacific 
profile  in  a parabolic  equation  model  (P.E.).  Tliose  calculations  pre- 
sented in  Fig.  III-2  have  also  been  computed  with  the  P.E.  model.  These 
results  are  presented  in  Fig.  III-5.  Comparisons  between  these  two  figures 
show  1 to  2 dB  differences  between  these  results  from  the  two  models. 
However,  the  results  from  both  models  do  have  similar  characteristics. 
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TABLE  III-l 


COMPARISON  BETWEEN  BOTTOM  isOUNCE  ENERGY  AND  BOTTOM  BOUNCE 
CONTRIBUTIONS  TO  THE  AVERAGE  PROPAGATION  LOSS 


Bottom  Loss 
Critical  Angle 
(deg) 


Bottom  Bounce  Energy 
(total  energy/water 
refracted  energy) 


Bottom  Bounce 
Contributions 
to  Prop.  Loss 


0 

0 

0 

1 

0.02 

0.25 

5 

0.47 

0.63 

10 

1.46 

1.50 

15 

2.50 

2.61 

20 

5.45 

5.47 
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TABLE  III-2 


COMPARISON  BETWEEN  BOTTOM  BOUNCE  ENERGY  AND  BOTTOM  BOUNCE 
CONTRXBOTIONS  TO  THE  AVERAGE  PROPAGATION  LOSS 


Source  Depth 
(m) 

Bottom  Bounce 

Energy 

Contributions  to 
Propagation  Loss 

53.5 

11.9 

12.2 

50.0 

5 . ‘r 

6.3 

90.0 

v.] 

3-6 

152.0 

• 3.’^ 

3.5 

600.0 

2.i4 

0.5 

-q8. 
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In  both  cases  the  water  column  is  divided  into  three  re^^ions.  'ilie  influences 
of  the  bottom  are  receiver  depth  dependent  in  the  first  and  third 
regions;  generally  the  sensitivity  to  bottom  loss  variation  Increases  as 
the  receiver  depth  approaches  the  surface  or  bottom.  Within  the  second 
region  the  bottom  influence  is  regular:  independent  of  the  receiver  depth. 
The  P.E.  model,  however,  predicts  additional  bottom  Influences  within  400 
to  700  m oT  the  bottom. 

Additional  P.E.  runs  were  made  using  various  bottom  depths.  The 
bottom  depths  used  were  364-7  m (critical  depth  minus  305  in),  3952  m 
(critical  depth),  4257  m (critical  depth  plus  305  n>),  and  5543  m (actual 
depth  for  this  profile).  Bottom  critical  angles  of  1°  and  10°  were  used 
to  determine  the  bottom  influence  in  a manner  similar  to  Figs.  III-3  and 
III-4.  (The  10°  angle  was  chosen  on  run  time  considerations  and  the  1° 
was  chosen  to  ensure  that  the  same  logic  was  used.)  The  differences  ob- 
tained from  the  results  using  these  two  bottom  descriptions  are  presented 
in  Fig.  III-6.  The  bottom  Influences  observed  from  Figs.  III-3  and  III-4 
are  again  present  in  Fig.  III-6.  An  additional  characteristic,  however, 
appears  in  the  very  near  bottom  region  (400  to  'JOO  m off  the  bottom).  In 
this  region  the  bottom  Influences  are  more  complex  and  are  also  dependent 
on  the  bottom  depth. 

1.  Additional  Profiles 


Profiles  from  the  North  Pacific  Ocean  and  from  the  Indian  Ocean 
were  also  used  in  this  study.  The  North  Pacific  profile  presents  a shallow 
sound  channel  axis  (120  m)  and  a 26  m/sec  sound  speed  difference  between 
the  surface  and  the  axis.  The  Indian  Ocean  profile  has  an  axis  depth  of 
1800  m and  a soimd  speed  difference  between  surface  and  axis  of  44  m/sec. 

The  North  Pacific  profile  was  used  to  study  bottom  influence  as  a function 
of  source  depth  (Fig.  III-7)  and  as  a function  of  bottom  depth  (Fig.  III-8) . 
The  Indian  Ocean  profile  was  used  to  present  the  bottom  influence  at 
various  bottom  depths  for  two  source  depths  (Figs.  III-9  and  III-IO). 
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The  source  depths  used  with  the  North  Pacific  profile  had 
conjugate  depths  of  19^5  m (50  m) , 151U  m (60  m),  850  m (76  m),  and  IO7  m 
(152  m) . In  Fig.  III-7  the  bottom  Influence  again  Increases  as  the  re- 
ceiver depth  approaches  either  the  surface  or  the  bottom  In  the  first  and 
third  regions  defined  by  the  soiarce  and  Its  conjugate  depth.  In  the  region 
between  these  two  depths,  the  bottom  Influence  no  longer  Is  as  well  de- 
fined as  It  is  for  the  mid-Paciflc  profile.  However,  the  propagation 
loss  within  this  region  attempts  to  present  a regular  dependence  upon  the 
bottom  loss. 

The  studies  on  the  effects  of  bottom  depth  variations  were 
conducted  using  the  P.E.  model  for  both  profiles.  For  the  North  Pacific 
profile  the  effects  on  propagation  loss  of  varying  bottom  depths  appears 
to  increase  the  region  of  a regular  dependence  upon  the  bottom  loss  down 
into  the  third  region  defined  by  the  source  conjugate  depth.  The  bottom 
influences  again  increase  as  a receiver  above  the  source  depth  approaches 
the  surface. 

The  two  source  depths  used  with  the  Indian  Ocean  profile  had 
conjugate  depths  of  3585  in  (152  m)  and  ^795  (18  m) . The  bottom  depths 
used  were  critical  depth  (^795  o')  > critical  depth  plus  and  minus  3^5  ™ 

(5100  m and  4490  m),  and  the  measiored  depth  (380O  m) . With  the  Indian 
Ocean  profile  the  bottom  influences  are  again  seen  to  be  divided  into 
three  regions.  The  major  effects  of  bottom  depth  variations  occur  within 
the  first  and  third  regions  along  with  the  very  near  bottom  effects  seen 
earlier.  However,  when  the  source  conjugate  depth  is  within  this  very 
near  bottom  region,  the  near  bottom  effects  extend  into  the  region 
between  the  source  and  the  source  conjugate  depth. 

2 . Bottom  Bounce  Energy  Only 

Further  investigations  have  included  modifications  to  FACT  to 
provide  propagation  loss  for  bottom  bounce  rays  only.  The  purpose  of 
these  modifications  was  to  determine  the  distribution  of  bottom  bounce 
intensity  throughout  the  water  column.  The  average  propagation  losses 
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were  calculated  and  plotted  in  the  same  manner  as  before.  Only  the  mid- 
Pacific  profile  has  been  used. 

Figure  III-ll  is  a representative  set  of  these  average  propagation 
loss  curves  for  the  bottom  bounce  rays.  The  bottom  loss  descriptions  for 
this  figure  are  again  a 0 dB  loss  below  the  critical  angle  and  a 50  dB 
loss  above.  The  distribution  of  bottom  boiince  intensities  (the  shape 
of  these  curves)  is  seen  to  be  virtually  independent  of  these  bottom 
descriptions.  (The  propagation  loss  for  the  1°  critical  angle  case  may 
have  been  undersampled  in  range,  resulting  in  the  modified  shape  of  the 
curve).  The  distribution  of  these  intensities,  althoxagh  not  presented 
here,  were  also  seen  to  be  indepiendent  of  variations  in  sovirce  depth. 

Figure  III-ll  could  have  represented  variations  in  source  depth.  Varia- 
f-lons  in  bottom  depth  simply  displaced  the  near-bottom  propagation  loss 
decay  in  depth,  leaving  the  shape  of  the  upper  portion  unchanged.  For 
these  bottom  loss  descriptions  the  distribution  of  bottom  bounce  energy 
appears  to  be  a function  of  the  profile  only.  The  bottom  bounce  energy 
is  within  5 dB  of  being  uniformly  distributed  throughout  the  water  column. 

Additional  bottom  loss  descriptions  were  also  investigated.  When 
the  bottom  loss  was  monotonically  increasing  with  grazing  angle  (Fig.  III- 
12),  the  propagation  loss  ctirves  again  were  characterized  by  the  shape 
of  the  curves  in  Fig.  III-ll.  Again,  variations  in  source  depth  and 
bottom  depth  did  not  greatly  modify  the  shape. 

The  bottom  loss  curve  (solid  line)  in  Fig.  III-IJ  was  also  used 
in  part  of  this  study.  Figure  III-14  presents  the  propagation  loss  curves 
for  four  source  depths.  Although  these  curves  are  still  similar  in  shape, 
there  are  definite  differences.  As  the  source  depth  increases,  the  near- 
axis  propagation  loss  increases. 

Two  modifications  to  this  bottom  loss  curve  were  used  (dashed 
and  dotted  curves  in  Fig.  III-lj)  to  isolate  the  cause  of  the  differences 
mentioned  above.  The  propagation  loss  curves  for  the  600  m source  appear 
in  Fig.  Ill- 15.  The  modified  bottom  loss  cvirve  represented  by  the  dotted 
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curve  generated  a curve  with  the  familiar  shape  presented  in  Fig.  III-ll. 
The  modified  bottom  loss  curve  represented  by  the  dashed  line  magnifies 
the  variation  in  the  shape  of  the  propagation  loss  curves  seen  in  Fig. 
III-14. 


The  propagation  loss  versus  range  ctirves  for  a monotonically 
increasing  bottom  loss  curve  and  for  the  bottom  loss  curve  in  Fig.  III- 
15  are  presented  in  Fig.  III-I6  for  the  6OO  m to  600  m source-receiver 
depth  combination.  This  comparison  was  made  to  help  isolate  the  cause  of 
the  differences  seen  in  Figs.  III-14  and  III-I5.  The  monotonically  in- 
creasing bottom  loss  results  in  a smooth  propagation  loss  curve.  The 
notched  bottom  loss  ciirve,  however,  results  in  a sawtooth  curve;  this 
indicates  there  is  a bundling  of  bottom  bounce  energy.  As  the  receiver 
moves  off  axis,  the  sawtooth  structure  decreases  in  amplitude  approaching 
the  smooth  structure  of  the  monotonically  increasing  bottom  loss. 

C.  Conclusions 

Bottom  influences  on  propagation  loss  divide  the  water  column  into 
three  regions  based  on  the  source  and  so\irce  conjugate  depths.  The  total 
propagation  loss  within  the  region  between  the  source  and  its  conj\igate 
depth  has  a regular  dependence  on  the  bottom  loss.  In  the  remaining  two 
regions  the  bottom  influence  is  dependent  on  the  receiver  depth.  As  the 
receiver  approaches  the  surface,  as  well  as  the  bottom,  the  bottom  in- 
fluence generally  increases.  In  the  very  near  bottom  region  (400  to  TOO  m 
from  the  bottom),  wave  theory  may  be  required  to  adequately  predict  the 
depth  dependence  of  the  bottom  Influence. 

Most  bottom  loss  curves  used  in  ray  theory  models  are  described  by 
a monotonically  increasing  reflection  coefficient.  Additional  attention 
needs  to  be  placed  on  anomalous  dips  seen  at  the  lower  grazing  angles  in 
bottom  bounce  data. 
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FIGURE  III-16 

PROPAGATION  LOSS  versus  RANGE  FOR  BOTTOM  BOUNCE  ENERGY  ONLY 
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IV.  EFFECTS  OF  RANGE  VARIABLE  BATHYMETRY 


This  chapter  contains  a synopsis  of  work  carried  out  during  this 
contract  year  on  the  problems  of  sloping  bottom  and  rough  interface  effects 
in  propagation. 

A.  Sloping  Bottom  Effects 

In  the  remainder  of  this  section  we  will  examine  some  data  which 
show  slope  enhancement.  The  data  will  be  compared  to  transmission  loss 
curves  generated  using  the  parabolic  equation  model. 

We  shall  define  slope  enhancement  as  the  increase  in  sound  level  over 
the  sloping  area  of  a reflecting  bottom  above  the  so\ind  level  measured 
over  the  same  sloping  area  of  an  absorbing  bottom  (critical  angled”). 

This  increase  is  Illustrated  in  Figs.  IV-l  and  IV-2.  The  transmission 
loss  curves  presented  in  these  figures  were  computed  using  the  parabolic 
equation  model  described  in  a previous  report  (Hawker  et  al.^).  The  bathy- 
metry and  sound  speed  profile  used  in  these  model  runs  are  shown  in  Fig. 
IV-3.  Computations  were  done  at  a frequency  of  93  Hz. 

During  a Naval  exercise  in  the  Atlantic  Ocean,  a continuous  wave 
source  was  driven  through  deep  water,  over  the  continental  slope,  and  over 
the  continental  shelf.  Figiire  IV-4  sketches  the  bathymetry  of  the  exer- 
cise area.  Note  that  the  steepest  part  of  the  continental  slope  lies 
between  the  ranges  of  175  and  I90  nml. 


Figure  IV-5  shows  a typical  sound  si)eed  profile.  All  of  the  sound 
speed  profiles  in  the  area  of  interest  were  very  similar  to  Fig.  IV-5. 

The  receiver  depths  having  good  data  are  indicated  in  the  main  section  of 
the  figure  and  the  source  depths  in  the  inset.  There  is  one  receiver  in 
each  sound  channel,  one  just  below  critical  depth,  and  one  on  the  upper 
border  between  the  two  sound  channels.  The  shallow  source  lies  in 
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the  surface  duct  and  the  deeper  source  lies  in  the  upper  sound  speed 
channel.  Transmission  loss  was  measured  at  35  and  93  Hz  for  the  shallow 
source  and  at  IO5  and  l44  Hz  for  the  deeper  source. 

The  data  have  been  smoothed  by  taking  an  equally  weighted  average 
of  observations  within  5 nmi  of  the  desired  range  point.  There  is  one 
bad  data  point  at  about  l4o  nmi  for  the  receiver  at  181O  m;  this  is 
suspected  to  be  an  unremoved  shot.  Such  a point  causes  a spike  in  the 
smoothed  data  spread  over  10  nmi;  this  area  is  shaded  out  of  the  data 
presented  in  Figs.  IV-6,  IV-7>  and  IV-8. 


As  a reference  device,  the  dotted  lines  on  Figs.  IV-6,  IV-7,  and  IV-8 
indicate  transmission  loss  given  by  simple  cylindrical  spreading.  These 
reference  lines  were  set  by  selecting  a source  level  that  would  match  the 
roughly  averaged  data  aroiind  100  nmi  and  computing  cylindrical  spreading 
from  this  selected  pseudo-original  source  level.  Note  that,  in  Figs. 

IV-1  and  rV-2,  the  transmission  loss  for  the  absorbing  bottom  drops 
sharply  after  the  leading  edge  of  the  continental  slope  is  encountered. 
Therefore,  data  levels  remaining  nearly  at  or  above  the  cylindrical 
spreading  line  throughout  the  slope  area  will  be  interpreted  as  indicating 
slope  enhancement. 

Using  the  criterion  discussed  above,  all  of  the  data  presented  show 
some  slope  enhancement.  The  data  at  95  Hz  (Fig.  rv-6)  is  clearest  and 
has  the  most  obvious  slope  enhancement  of  all  the  data  examined.  Trans- 
mission loss  is  enhanced  by  approximately  8 dB  above  the  cylindrical 
spreading  line  for  all  but  the  deepest  hydrophone.  The  deepest  (2467  m) 
hydrophone  is  below  critical  depth  and  evidences  reduced  slope  enhancement 
as  well  as  a generally  lower  level  for  all  frequencies  and  source  depths 
examined  here. 

The  data  taken  at  IO5  Hz  (Fig.  IV-7)  shows  a moderate  slope  enhancement 
of  about  5 dB  above  the  cylindrical  spreading  line  for  the  715 
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and  1810  m hydrofiiones.  The  hydrophone  at  5^5  shows  less  enhancement 
about  5 than  the  middle  hydrophones. 

At  l44  Hz  (Fig.  IV-8)  the  585  ra  hydrophone  shows  moderate  slope 
enhancement  and  1810  m hydrophone  strong  enhancement.  The  715  hydrophone 
data  evidence  less  enhancement  than  the  data  of  other  frequencies  at  this 
depth.  The  2467  m hydrophone  at  l44  Hz  shows  less  slope  enhancement  than 
any  of  the  data.  The  data  average  to  the  cylindrical  spreading  level  but, 
by  the  criterion  stated  earlier,  slope  enhancement  is  still  present. 

The  35  Hz  data  reflected  cable  strumming  contamination  and  is  not 
considered  in  this  report. 

In  a previous  report,  Hawker  et  al,^  drew  three  conclusions  from 
the  modeled  transmission  loss  curves  then  available.  We  will  compare  the 
conclusions  quoted  from  the  report  with  the  data  presented  here. 

(1)  Model  conclusion  — All  receiver  depths  show  a slope  enhance- 
ment with  enhancement  increasing  with  decreasing  depth. 

Data  conclusion  --  For  the  modeled  transmission  loss  curves, 
"source  depth"  is  receiver  depth  for  the  data  and  vice  versa.  Rectifying 
the  switch  in  terms,  shallower  sources  show  more  slope  enhancement  for 
the  modeled  curves.  All  data  receiver  depths  show  slope  enhancement.  The 
data  from  the  9I  source  shows  uniformly  less  enhancement  than  the  data 
from  ^he  18  m source.  The  change  in  frequency  from  93  Hz  to  105  Hz  should 
not,  by  itself,  cause  any  significant  differences. 

(2)  Model  conclusion  --  The  enhancement  begins  at  a location 
corresponding  to  approximately  40^  up  the  continental  slope  and  peaks  at 
the  top  of  the  slope. 

Data  conclusion  --  The  slope  enhancement  becomes  apparent  at 
about  4o^  up  the  continental  slope,  a range  of  about  I80  nmi.  The  enhancement 
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peaks  at  or  very  near  the  top  of  the  slope  in  all  cases. 


(3)  Model  conclusion  --  Although  the  enhanced  energy  decays  with 
distance  beyond  the  top  of  the  slope,  the  decay  rate  is  nearly  the  same 
in  all  cases. 


Data  conclusion  --  For  all  data  at  one  frequency,  the  slope 
enhancement  decays  at  nearly  the  same  rate  to  nearly  the  same  level  for 
all  receiver  depths.  The  deepest  hydrophone  decays  at  the  same  rate  as 
its  shallower  compatriots  and  to  a level  decreased  by  4 to  8 dB. 

The  data  conclusions  generally  agree  with  those  drawn  from  the  modeled 
transmission  loss  curves.  The  models  show  promise  as  an  effective  tool 
for  studying  slope  enhancement  aspects  of  the  sloping  bottom  problem. 

In  addition  to  the  work  just  described,  additional  studies  have  been 
carried  out  on  the  feasibility  of  employing  normal  mode  coupling  methods 
to  the  solution  of  sloping  bottom  problems.  The  chief  advantage  of  such 
techniques  would  be  the  automatic  incorporation  of  a detailed  description 
of  the  bottom  into  the  calculations.  Assessments  of  existing  techniques, 
such  as  ray  trace  methods  and  a parabolic  equation  approach  already  carried 
out  under  this  program,  indicate  that  it  is  necessary  to  consider  the  sub- 
bottom structure  in  detail,  especially  at  very  low  frequencies. 

The  work  carried  out  on  a mode  coupling  approach  is  still  in  a 
formative  stage  and  no  quantitative  results  are  yet  available.  It  is 
expected  that  during  the  next  contract  year  (FY  7?)  a definite  decision 
will  be  made  concerning  the  applicability  and  usefulness  of  this  approach. 

B.  Rough  Interface  Studies 

The  general  goal  of  the  rough  interface  studies  being  conducted  as 
part  of  the  Naval  Electronic  Systems  Command  (NAVELEX),  Code  320,  sponsored 
bottom  Interaction  program  is  to  assess  the  usefulness  of  existing 
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techniques  in  accounting  for  the  influence  of  bottom  roughness  on  propagation. 
Also,  the  most  viable  of  these  techniques  are  to  be  used  in  a sensitivity 
study  designed  to  define  roughness  regimes  which  are  of  practical  importance 
to  propagation  problems. 

During  the  present  contract  year  (fY  76),  research  efforts  in  the 
area  of  rough  interface  effects  have  centered  around  the  theoretical 
understanding  of  scattering  theory  and  its  applicability  to  propagation 
problems,  as  well  as  the  natui’e  and  usefulness  of  a smoothed  (or  average) 
boundary  condition  approach.  These  efforts  are  still  largely  in  a forma- 
tive stage  and  will  not  yie].d  significant  quantitative  results  until  the 
next  contract  year  (FY  77)- 

Tlie  remainder  of  this  section  is  devoted  to  a discussion  of  the  nature 
of  the  rough  surface  problem  and  its  relationship  to  the  general  problem 
of  range  changing  bathymetry. 

A plane  wave  incident  on  a statistically  rough  surface  at  a given 
angle  gives  rise  to  an  outgoing  scattered  wave  which  is  made  up  of  an 
average  or  coherent  field  and  a fluctuating  or  incoherent  field.  If  the 
surface  is  only  slightly  rough  (on  the  scale  of  a sound  wavelength),  the 
scattered  energy  will  be  primarily  coherent  and  will  be  concentrated  in 
the  specular  direction  (reflection  from  a lossy  mirror).  As  the  roughness 
increases,  the  coherent  component  decreases  and  the  incoherent  :omponent 
dominates,  causing  the  scattered  field  to  become  more  and  more  diffuse.  In 
a propagation  context  these  processes  would  be  described  in  terms  of  mode 
conversion  and  energy  flow  from  the  coherent  field  into  the  incoherent 
field.  Mathematically,  the  two  problems,  scattering  and  propagation,  are 
quite  different.  In  the  scattering  problem  one  seeks  an  asymptotic  quantity, 
a scattering  coefficient;  in  a propagation  problem,  however,  the  roughness 
affects  the  mode  structure  everywhere  and  nearfleld  effects,  as  well  as 
coupling  to  refractive  effects,  cannot  be  ignored. 

It  is  clear  from  research  on  scattering  theory  that  both  the  amplitude 
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and  phase  of  scattered  fields  can  be  appreciably  altered  by  small  roughnesses. 
This  will  cause  a modification  of  arrival  angle  structure  seen  by  an  array  as 
well  as  a redistribution  of  energy  along  the  propagation  path  in  such  a 
way  that  the  propagation  loss  is  increased. 

Various  geophysical  processes  shaping  the  ocean  basins  have  led  to 
a wide  spectrum  of  bottom  roughness  scales,  ranging  from  the  walls  of  the 
ocean  basins  (continental  slopes)  to  minute  ripples  on  otherwise  smooth 
abyssal  plains.  Figure  IV-9  shows  a schematic  representation  of  these 
scales  as  well  as  the  ranges  over  which  echo  sounding  and  other  methods 
can  be  used. 

The  bottom  interaction  program  at  ARL:UT  has  concentrated  on  examining 

the  effects  of  bathymetry  variations  having  wave  numbers  on  the  extremes 

of  this  spectrum.  On  the  one  hand  there  are  the  very  large  scale  effects 

10^  m to  10^  m which  are  being  investigated  under  the  heading  of  sloping 

bottom  problems.  Examples  of  effects  in  this  range  were  given  in  the 

previous  two  sections  of  this  chapter.  At  the  other  extreme  there  are  the 

small  scale  irregularities,  treated  stochastically,  which  are  part  of  the 

bottom  roughness  study.  As  ever,  the  intermediate  ground  is  the  most 

difficult  and  techniques  applied  to  either  extreme  are  not  necessarily 

well  adapted  for  use  there.  The  topic  of  roughness  scales  as  well  as 

2 

acoustic  interaction  is  discussed  at  length  by  Clay  and  Leong. 

The  points  made  about  topograjiilcal  scales  in  Fig.  IV-9  can  be  given 
added  force  with  a realistic  example.  The  bathymetry  of  the  baseline 
track  for  a Naval  exercise  in  the  Pacific  Ocean  was  analyzed  for  its 
statistical  properties.  The  distribution  of  heights  around  the  mean  plane 
for  this  section  is  given  in  Fig.  IV-10.  The  rms  height  was  found  to  be 
119  “>•  The  correlation  coefficient  of  heights  given  in  Fig.  IV-11  is  seen 
to  decay  rapidly  toward  zero  at  approximately  25  km,  followed  by  quasi- 
perlodic  behavior.  The  rms  slope  computed  for  this  track  interval  is 
1.1’’.  The  "average  length"  computed  using  an  "average"  slope  of  1.1°  and 
an  "average"  height  of  120  m is  6.25  km,  which  agrees  well  with  the  range 


-130- 


RELIEF 


1 

i GROSS  FEATURES 

INTERMEDIATE  FEATURES 

SMALL  FEATURES 

MID  OCEAN  RIDGES 

BANKS 

RIPPLES 

ABYSSAL  PLAINS 

VALLEYS 

BOULDERS 

1 TRENCHES 

HILLS 

MOUNDS 

LARGE  HILLS,  etc. 

etc. 

ROCK  OUTCROPS 

etc. 

FIGURE  32:-  9 

SPECTRUM  RANGE  OF  SEA  FLOOR  ROUGHNESS  AND  RESOLUTION 
CAPABILITIES  OF  MEASUREMENT  TECHNIQUES 


fAdopted  from  Ako 


25, 


131 


ARL - UT 
AS-77-279 
KEH  - DR 
3 -22  -77 


TRUE  PROBABILITY  0.073 


132 


HEIGHT  ABOVE  MEAN  DEPTH  (m) 

FIGURE  IV-10 

RELATIVE  PROBABILITY  DENSITY  OF  HEIGHTS  ABOVE  MEAN  DEPTH  FOR  CHURCH  ANCHOR 
BASELINE  TRACK  FROM  SITE  B TO  SITE  C 


CORRELATION  COEFFICIENT 


o 


Range  - km 
FIGURE  IV-n 


CORRELATION  COfePFICIENT  OF  WATER  DEPTHS  versus  RANGE 
ALONG  CHURCH  ANCHOR  BASELINE  FROM  SITE  B 


133 


AS-77-288 


at  which  the  correlation  function  shown  in  Fig.  IV-2  has  decayed  to  a value 
of  e (one  definition  of  a correlation  length).  The  resulting  wave 
number,  2m/6.25  ton  = ;lo  ^ m and  height,  120  m,  leads  to  an  "average 
feature"  which  lies  nearly  on  the  solid  line  of  Fig.  IV-9>  near  the 
boiuidary  separating  large  features  from  intermediate  features.  Thus, 
this  abyssal  hills  region  is  one  which  lies  just  at  the  outer  bovinds  of  the 
sloping  bottom  problem. 

There  are  a variety  of  acoustic  and  nonacoustic  techniques  which  can 
be  used  to  obtain  information  about  bathymetry.  The  shaded  areas  in  Fig. 
IV-9  show  the  domains  in  which  conventional  vertical  echo  sounders,  deep 
towed  high  resolution  systems,  and  underwater  photography  are  operative. 
Conventional  vertical  echo  sounders,  having  beamwidths  of  approximately 
10°,  span  the  large  and  intermediate  feature  regions  of  the  spectrum 
adequately.  However,  the  small  feature  region  where  the  rou^ness  problem 
lies  is  Inaccessible  by  these  conventional  techniques.  Underwater  photo- 
graphy is  both  time  consuming  and  expiensive  and  clearly  not  well  adapted 
to  wide  area  coverage.  It  would  seem,  therefore,  that,  if  roxighness  is 
to  be  taken  into  account,  the  bottom  roughness  information  will  have  to 
come  from  some  other  source,  such  as  acoustic  experiments. 

The  bottom  roughness  problem  is  thus  one  for  which  one  gets  no 
intuitive  feeling  by  examining  ordinary  bathymetry  charts.  The  goal  of 
research  in  this  area,  particularly  the  sensitivity  studies  being  conduct- 
ed at  ARL:UT,  is  to  develop  an  approximate  method  for  accounting  for  the 
acoustical  effects  of  such  small  features  without  requiring  the  use  of 
detailed  statistical  Information  which  may  not  be  available. 

This  effort  is  based  on  the  use  of  a smoothed  boundary  condition 
applied  to  the  mean  (flat)  surface.  Initial  work  suggests  that  this 
approach  can  be  turned  into  a well-defined,  useful  method  for  obtaining 
roughness  effects  in  which  one  can  have  high  confidence. 

Although  it  is  by  no  means  clear  that  any  results  from  rough  surface 
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scattering  theory  can  be  used  to  provide  boundary  conditions  in  a propagation 
model,  it  is  nevertheless  useful  to  examine  one  case  of  rough  surface 
scattering.  Under  reasonably  general  circumstances  it  can  be  shown  that 
the  intensity  of  the  coherent  scattered  field,  I^,  can  be  given  in  terms 
of  the  field  reflected  from  a plane  stirface  of  the  same  material,  I^,  by 
the  expression  I = a , I , where  a , is  the  coherent  scattering  coefficient. 
Fig\ire  IV-12  shows  an  example  of  the  comparison  of  computed  coherent 
scattering  coefficients  with  measurements  made  at  ARL:UT.  These  data  have 
been  reported  by  Boyd  and  Deavenport,^  The  Rayleigh  parameter  g is  defined 
as  4k  h sin  6,  where  k is  the  acoustic  wave  number,  h the  rms  height,  and 
6 the  grazing  angle.  It  will  be  observed  that,  for  small  roughness  (small 
g),  both  theoretical  predictions  converge  on  the  data  and  agreement  is 
good.  For  the  range  0 < g < 1.5  the  scattering  can  be  expected  to  be 
primarily  coherent  and  theoretical  prediction  of  the  scattering  coefficient 
will  be  good.  It  is  just  this  regime  which  is  important  in  the  low  (and 
very  low)  frequency  bottom  interaction  problem.  For  example,  at  100  Hz 
a feature  1 m high  will  produce  a Rayleigh  parameter  value  less  than  O.7O. 

Such  a featiire  with  a mean  length  of  6 m will  lie  near  the  center  of  small 
roughness  regime  of  Fig.  IV- 9. 
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FIGURE  I3Z:-12 
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V.  ANALYTICAL  TECHNIQUES  FOR  DETERMINING  SUBBOTTOM 
VELOCITY  PROFILES  IN  UNCONSOLIDATED  SEDIMENTS 


The  purpose  of  this  chapter  is  to  discuss  and  simmarize  research 
previously  released  in  an  ARL:UT  technical  report  by  Rutherford.^  The 
purpose  of  this  work  was  to  describe  and  summarize  different  analytical 
methods  for  determining  velocity  profiles  in  unconsolidated  sediments  from 
experimental  data.  This  work  was  by  no  means  exhaustive;  its  intent  was 
to  assemble  in  one  place  those  methods  which  would  be  most  useful  to  under- 
water sound  propagation  modelers  and  best  suited  to  the  types  of  data  to 
which  they  usually  have  access. 


There  are  different  ways  of  obtaining  information  about  the  velocity 

structure  of  the  subbottom.  In  shallow  water,  divers  with  probes  can  take 

samples  of  the  upper  few  meters  of  sediment  for  later  laboratory  analysis. 

A somewhat  deeper  probing  of  the  sediment  structure  is  possible  with 

sediment  coring  apparatus  deployed  from  a research  vessel.  The  cores  thus 

obtained  may  be  analyzed  in  the  laboratory  or  the  velocity  profile  of  the 

sediment  may  be  measured  in  situ  with  a device  such  as  the  ARLrUT  Compres- 

2 

sional  Wave  Profilometer , which  can  be  attached  to  a corer.  These  direct 
methods  are  limited  by  the  depths  that  a corer  can  probe,  roughly  ^0  m. 

The  only  other  direct  methods  involve  the  analysis,  either  in  situ  or 
laboratory,  of  drilling  cores.  This  method,  which  is  very  useful  in  many 
acoustic  studies,  relies  on  well  cores  which  are  often  widely  spaced  and 
must  be  supplemented  with  other  methods. 


The  methods  addressed  in  the  abovementioned  report  are  indirect, 
acoustical  techniques,  which  more  often  than  not  provide  the  only  way  of 
estimating  the  velocity  structures  of  the  unconsolidated  sediments  to  the 
depths  required  by  underwater  propagation  modelers.  These  indirect  methods 
are  concerned  only  with  the  unconsolidated  sediments,  although  some  of  the 
methods  can  be  and  are  used  in  seismic  profiling.  Seismological  methods 
such  as  those  Involving  lateral  and  shear  waves  were  not  discussed  in  the 
report.  The  report  was  primarily  Intended  to  serve  the  needs  of  those 
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working  in  the  field  of  underwater  sound  propagation  modeling  who  have 
little  knowledge  of  the  techniques  for  extracting  subbottom  velocity  in- 
formation which  is  necessary  input  for  their  models. 

The  remaining  sections  of  this  chapter  will  summarize  the  different 
methods  that  were  discussed  and  detailed  in  the  abovementioned  report. 

The  intent  of  these  methods  is  to  extract  information  concerning  sound 
speed  gradients,  layering  structure,  and  layer  velocities  from  experimental 
data,  primarily  wide  angle  profiling  data. 

The  sximmary  in  the  following  sections  of  this  chapter  deals  with  the 
different  methods  in  the  order  they  were  discussed  in  the  original  report. 
The  different  analytical  techniques  treated  are  the  reflected  arrival 
methods,  the  V(t)  method,  and  the  refracted  arrival  methods.  In  the 
original  report,  the  treatment  of  each  method  comprised  an  entire  chapter. 
In  the  following,  the  discussion  of  the  analytical  methods  will  be  drama- 
tically condensed.  Each  method  will  be  briefly  described,  the  results  of 
any  examples  will  be  quoted,  and  one  example  showing  the  connection  between 
the  reflected  arrival  methods. 'and  the  V(t)  method  will  be  Included. 

A.  Reflected  Arrival  Methods 

The  first  analytical  techniques  treated  sire  the  reflected  arrival 
methods.  These  methods  are  most  useful  when  wide  angle  subbottom  profiling 
data  which  exhibit  a sediment  layering  structure  are  analyzed.  The  pur- 
pose of  the  reflected  arrival  methods  is  to  determine  the  sound  velocities 
and  thickness  of  subbottom  sediment  layers  from  travel  time  curves  obtained 
from  wide  angle  profiling. 

There  are  two  basic  types  of  reflected  arrival  methods,  the  !I^(X) 

X 4 b 6 7 

methods  and  the  ray  parameter  methods.  In  the  derivation  of  both 

of  these  methods,  it  is  necessary  to  make  certain  simplifying  assumptions 

about  the  velocity  structure  of  the  layers  being  analyzed  so  that  tractable 

mathematical  results  may  be  obtained.  The  validity  of  these  simplifying 
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asslijnptions  ultimately  determines  how  well  the  T (X)  and  the  ray  parameter 
methods  perform. 

1.  The  T^(X)  Method^^^’^-^ 

The  three  T^(X)  methods  are  the  reduced  time-reduced  distance 
method,  the  short  range  method,  and  Durbaum's  method.  The  methods  are 
based  on  the  fact  that  the  relationship  between  the  travel  time  of  a ray 
reflecting  in  a single  isovelocity  layer  and  the  range  attained  while 
the  ray  is  traveling  through  that  layer  is  a hyperbolic  one  given  by 

T^(X)  = 1^(0)  + -^  . 

/-I  ^ 


The  quantity  denoted  by  is  the  isovelocity  speed  of  so\znd  in  the  layer 
and  T(0)  is  the  vertical  two-way  time  through  the  layer.  The  linear  re- 
lationship between  T^(X)  and  X^  for  the  isovelocity  layer  allows  and 
T(0)  to  be  obtained  graphically  if  the  travel-time  curve  (T  versus  X)  of 
arrivals  reflecting  in  the  layer  is  known. 

When  the  1^(X)  methods  are  applied  to  wide  angle  profiling  data, 
the  layers  exhibited  by  the  profiling  data  are  assumed  to  be  isovelocity. 
This  isovelocity  assumption  allows  one  to  analyze  travel-time  curves  of 
subbottom  layers  by  subtracting  the  effects  of  the  overlying  layers  such 
that  a one-layer  problem  remains.  When  this  is  done,  the  velocity  of 
sound  and  layer  thickness  in  a subbottom  layer  may  be  obtained  in  the  same 
manner  that  one  would  use  for  the  single  layer  case. 


The  isovelocity  assumption  necessary  in  the  derivation  of  the 
I^(X)  methods  has  two  major  effects.  First,  when  the  T^(X)  methods  are 
applied  to  nonlsoveloclty  situations,  the  sound  velocities  obtained  are 
interpreted  as  Interval  velocities,  which  are  the  time  averaged  sound 
speeds  through  the  layer.  Second,  the  methods  cannot  be  applied  to 

slt\iations  in  which  the  ratio  of  the  thickness  of  the  overlying  layers  to 
that  of  the  subbottora  layers  being  analyzed  is  greater  than  about  15:1. 
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This  problem  arises  because  the  absolute  error  introduced  by  the  isovelocity 
assumption  in  overlying  layers  is  magnified,  in  the  layer  being  analyzed, 
by  the  ratio  of  the  thickness  of  the  overlying  layers  to  the  thickness  of 
the  layer  in  question.  When  this  ratio  exceeds  about  15:1,  the  'I^(X) 
methods  can  no  longer  be  relied  upon  to  give  accurate  results. 


In  the  text  of  the  abovementioned  report  the  reduced  time-reduced 
distance  method  and  the  small  range  method  are  illustrated  with 
exeimples.  Both  of  these  methods  are  used  to  find  the  layer  velocity  in  a 
200  m thick  sediment  layer  underlying  a 500  m thick  water  layer  (thickness 
ratio  2.5:l).  The  reduced  time-reduced  distance  method  determined  a layer 
velocity  of  1559  m/sec  and  the  small  range  method  yielded  1575  m/sec.  The 
standard  value  calculated  from  equations  derived  in  the  report  was  I561 
ra/sec,  which  yielded  percentage  differences  of  l.^:-  and  O.89,  respectively. 

7 

2 . The  Ray  Parameter  Method 


The  ray  parameter  method,  which  is  the  second  of  the  reflected 
arrival  methods,  is  fundamentally  different  from  the  7^(X)  methods,  al- 
though it  also  analyzes  travel-time  data  of  reflected  arrivals.  The  ray 
parameter  method  is  unique  in  that  water  or  sediment  layers  overlying  the 
layer  to  be  analyzed  are  Irrelevant  to  the  solution;  thus,  simplifying 
assiamptions  about  the  velocity  structure  above  the  layer  in  question  need 
not  be  made.  The  only  simplifying  assumption  made  is  that  the  layer  being 
analyzed  is  isovelocity.  This  has  the  effect  of  allowing  the  ray  parameter 
method  to  be  used  in  layers  that  are  thinner  than  could  be  handled  by  the 
T^(X)  methods . 


The  ray  parameter  method  is  so  named  because  it  compares  the  ray 
paths  of  the  same  ray  parameter  which  have  reflected  off  the  top  and  bottom 
of  the  layer  in  question.  The  ray  parameter  is  defined  by  p = cos  9 (z)/c 
(z)  and  is  calculated  from  travel-time  data  by  evaluating  the  derivative 
of  a travel-time  curve,  l.e.,  p=dT/dx.  Two  different  types  of  ray  para- 
meter methods  are  discussed  in  the  report,  the  normal  method  and  the  thin 
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layer  approximation.  Both  methods  calculate  the  sound  velocity  and 
thickness  of  a particular  layer  by  matching  the  derivatives  of  travel- 
time curves.  The  thin  layer  approximation  incorporates  a graphical 
technique  which  alleviates  instabilities  associated  with  thin  layers  so 
that  thickness  ratios  of  up  to  250:1  may  be  handled. 

An  example  of  the  ray  parameter  method  is  considered  in  the 
text  of  the  above mentioned  report.  The  results  will  be  summarized  here. 

The  thin  layer  approximation  of  the  ray  parameter  method  was  used  on  a 
50  m subbottom  layer  having  a constant  sound  speed  gradient  underlying  a 
5000  m thick  layer  (thickness  ratio  100:1).  The  thin  layer  approximation 
calculated  a layer  thickness  of  49  m and  an  interval  velocity  of  1549 
m/sec.  This  Interval  velocity  compared  favorably  with  the  value  calculated 
from  an  equation  derived  for  a layer  having  a linear  velocity  profile 

( 1565  m/ sec) . 

B.  The  V(t)  Method^’ ^ 

The  next  technique  discussed  is  a method  whereby  velocity  measurements 
obtained  by  reflected  arrival  methods  may  be  analyzed  to  construct  a ve- 
locity profile  throughout  the  sediment  layers.  As  stated  earlier,  the 
reflected  arrival  methods  determine  a layer  velocity  under  the  isovelocity 
assumption,  a condition  which  is  rarely  met  in  the  physical  world.  For 
this  reason,  the  velocities  calculated  by  the  reflected  arrival  methods 
are  interpreted  as  interval  velocities  which  are  assmed  to  be  the  time 
averaged  sound  velocity  in  a particular  layer  or  interval. 

The  V(t)  method  is  a method  which  analyzes  the  depth  variation  of 
interval  velocities  to  obtain  velocity  gradients  and  layer  thicknesses 
throughout  the  sediment  layer  structiore.  The  method  is  based  on  a velocity 
function  which  is  formed  by  fitting  the  layer  or  interval  velocities  to 
a polynomial  in  t,  the  one-way  vertical  travel  time  to  the  time  midpoint 
of  the  various  layers  in  the  sediment  structure.  Once  the  velocity  function 
is  known,  a gradient  function  and  a depth  function  which  specify  the  gradient 
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and  depth  as  a function  of  t may  easily  be  determined. 

The  V(t)  method  is  implemented  by  fitting  the  interval  velocities  to 
a continuous  third  order  polynomial  in  t,  the  one-way  travel  time  to  the 
layer  midpoints,  i.e.. 


V(t)  = A + Bt  + Ct  + Dt' 


The  depth  as  a function  of  t (depth  function)  is  obtained  by  integrating 
V(t)  to  obtain 


/ 


H(t)  =#  V(t)dt  = At  + 


iBt^ 


i Ct^ 


1 


Dt 


The  velocity  function  may  be  manipulated  in  the  following  way  to  give  a 
gradient  function,  g(t). 


dV(t)/dt  dV(t)/dt 
~ dH(t)/dt  “ V(t) 


g(t) 


B + 2CT  + jPt 


A + Bt.  + Ct.  + Dt 


5 


The  gradient  and  depth  functions  may  then  be  combined  to  yield  the  gradient 
as  a function  of  depth. 

The  main  limitations  of  the  V(t)  method  are  associated  with  errors  in 
the  interval  velocities  and  velocity  discontinuities.  In  the  text  of  the 
report  which  this  chapter  summarizes,  these  effects  are  examined  and  shown 
to  be  significant  in  some  situations.  In  particular,  an  example  Involving 
velocity  discontinuities  predicted  layer  gradients  that  were  in  error  by 
an  average  of  50^  ss  opposed  to  average  differences  of  when  velocity 
discontinuities  were  absent. 

As  an  example  relating  the  reflected  arrival  methods  and  the  V(t)  method. 
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consider  the  case  of  four  subbottoni  layers  underlying  a 5000  m water  layer 
having  a bilinear  profile.  Table  V-1  lists  the  parameters  of  the  sub- 
bottom layers  and  the  water  layer.  The  subbottom  profile  of  this  example 
has  gradient  discontinuities  across  the  layer  Interfaces  but  the  sound 
velocities  are  continuous. 

Figure  V-1  depicts  the  computer  generated  T-X  cm’ves  for  the  example. 
Since  the  thickness  ratios  are  large,  the  ray  parameter  method,  using  the 
thin  layer  appT'oxlmation,  was  employed  to  determine  the  Interval  velocities, 
layer  thicknesses,  and  vertical  travel  times.  These  values  are  listed 
In  Table  V-2.  For  comparison,  the  actual  layer  thicknesses  and  the  calcu- 
lated layer  velocities  are  Included  In  parentheses. 

Before  applying  the  V(t)  method,  one  needs  to  form  a set  of  values 
for  the  Interval  velocity  versus  one-way  travel  time  to  the  midpoints.  In 
time,  of  the  four  subbottom  layers.  These  values  are  listed  In  Table  V-5. 
The  one-way  times  are  referenced  to  the  Interface  between  the  water  and  the 
first  subbottom  layer,  l.e.,  t=0  refers  to  the  top  of  the  first  subbottom 
layer . 

The  velocity-time  values  are  next  used  as  Input  to  the  polynomial 
fitting  procedure.  In  this  example,  the  values  of  Instantaneous  velocity 
at  the  top  of  the  first  layer,  l.e.,  t=0,  are  assumed  to  be  known  and  are 
Included  In  the  fitting  procedure.  One  output  of  this  procedure  Is  the 
velocity  as  a function  of  depth.  A plot  of  this  function  Is  pictured  In 
Fig.  V-2.  The  vertical  bars  In  Fig.  V-2  denote  typical  uncertainties  of 
±100  m/sec  In  the  Interval  velocities.  From  the  values  of  velocity,  depth, 
and  gradient  as  a function  of  t.  It  Is  possible  to  reconstruct  the  sub- 
bottom profile.  The  velocity  profile  determined  In  this  manner  for  the 
present  example  Is  listed  In  Table  V-4.  Again,  the  correct  values  are 
listed  in  parenthesis. 


TABLE  V-I 


WATER  PROFILE 


Depth 

Sound  Velocity 

(ra) 

(m/ sec) 

0 

1500 

1280 

1487 

5000 

1550 

SUBBOTTOM  PROFIIi: 


Subbottom  Thickness  Constant  Gradient  Initial  Sound 

Layer  No.  (m)  Velocity 

(m/si 


1550 

1600 


ARl.  . UT 
AS- 76-  1182 
SRR  - DR 
11-2-76 


FIGURE  I-l 

TRAVEL  TIME  CURVES  FOR  A THICK  WATER  LAYER 
OVERLYING  FOUR  THIN  SEDIMENT  LAYERS 


TABLE  V-2 


LAYER  PARAMETERS  CALCULATED  BY  THE  THIN  LAYER 
APPROXIMATION  OF  THE  RAY  PARAMETER  METHOD 


Subbottorn 
Layer  No. 

^Z 

(m) 

Interval  Velocity 
(m/ sec) 

2 -Way  Vertical  Travel 
Time 
(sec) 

1 

99  (100) 

1560  (1575) 

0.1270 

2 

101  (100) 

1660  (1650) 

0.1212 

5 

100  (100) 

1779  (177^) 

0.1127 

4 

101  (100) 

1957  (19^9) 

0.1027 
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TABLE  V-3 


INTERVAL  VEIDCITY  AS  A FUNCTION  OF  ONE-WAY  TIME 


Interval  Velocity 
(ra/ sec) 


One-Way  Time  to  Midpoint 
of  Layer 

( sec) 


1560 

0.03175 

1660 

0.09381 

1779 

0.15228 

1957 

0.20613 
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FIGURE  -SL-2 

SUBBOTTOM  VELOCITY  PROFILE  OBTAINED  BY  V(t)  METHOD 


ARL  - UT 
AS-76-n84 
SRR  - DR 
11-2-76 
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TABLE  y-k 


VELOCITY  PROFILE  RECONSTRUCTED  FROM  V(t)  METHOD 


Subbottom 
Layer  No. 

Thickness 

(m) 

Gradient 
(sec  ^) 

Initial  Sound 

Velocity 

(m/sec) 

1 

99.7  (100) 

0.5  (0.5) 

1550  (1550) 

P 

95.2  (100) 

1.1  (1.0) 

1602  (1600) 

5 

77.0  (100) 

1.5  (1.5) 

1716  (1700) 

4 

129.0  (100) 

1.8  (2.0) 

1862  (1850) 
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C.  The  Refracted  Arrival  Methods 


The  last  methods  developed  are  the  refracted  arrival  methods,  Dicus' 
method^^  and  Hanna's  method. These  methods  differ  from  the  reflected 
arrival  methods  in  that  the  effects  of  shallow  refracted  arrival  methods 
in  that  the  effects  of  shallow  refracted  arrivals  Instead  of  reflected 
arrivals  are  analyzed;  also,  these  methods  explicitly  assiame  that  a gra- 
dient exists  in  the  sediment  layer.  It  is  expected  that  both  of  these 
methods  would  most  successfully  be  applied  to  thick,  single  sediment  layers. 
The  term  refracted  arrival  should  not  be  confused  with  arrivals  resulting 
from  lateral  or  Interface  waves;  it  means  those  arrivals  that  are  returned 
from  the  sediment  layer  by  virtue  of  a turning  point  arising  because  of  a 
velocity  gradient. 

Dicus'  method  is  a graphical  method  which  deals  with  the  time  differences 
of  the  refracued  arrivals  and  the  arrivals  that  attain  the  same  receiver 
position  via  one  bottom  reflection.  The  method  is  derived  assuming  a 
pseudciinear  sediment  velocity  profile  and  yields  the  gradient  at  the  top 
of  the  sediment  layer  and  the  ratio  of  velocities  at  the  water  sediment 
interface  by  a graphical  solution  of  the  following  equation 


The  quantity  At  is  the  difference  in  travel  time  of  the  reflected  and 

refracted  arrivals.  C /C  is  the  ratio  of  the  sound  velocities  at  the 

o'  w 

water- sediment  Interface,  9^  is  the  bottom  angle  of  the  reflected  arrival, 
and  p is  the  gradient  at  the  top  of  the  sediment  layer. 

In  an  example,  Dicus'  method  was  applied  to  a single  sediment  layer 
with  a constant  gradient  of  1 sec~^  having  a velocity  ratio  at  the  water- 
sediment  Interface  of  1.0353 . Dicus'  method  calculated  a velocity  ratio 
of  1.0358  and  a gradient  at  the  top  of  the  layer  of  0.99  sec  which 
correspond  to  differences  of  0.24')t  and  1^,  respectively. 
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Finally,  a completely  different  method  by  Hanna  is  considered.  This 
technique  yields  the  same  output  as  Dicus'  method  by  an  analysis  of 
the  nulls  in  propagation  loss  curves  which  arise  because  of  the  interference 
produced  by  refracted  arrivals.  This  method  also  models  the  sediment  with 
a pseudolinear  velocity  profile. 

Hanna's  technique  proposes  to  determine  sediment  velocity  and  gradient 
information  from  the  interference  structure  exhibited  by  the  transmission 
loss  curve  for  a bottomed  receiver  over  the  range  of  the  direct  arrival. 

The  range  is  restricted  in  this  manner  because  outside  the  range  of  the 
direct  arrival,  multibounce  ray  paths  fiorther  complicate  the  transmission 
loss  by  adding  more  Interference  structure. 

If  there  were  no  energy  being  returned  from  the  bottom,  the  interference 
structure  of  the  transmission  loss  should  be  accounted  for  quite  well  by 
the  Lloyd's  Mirror  effect.  (The  Lloyd's  Mirror  effect  is  most  noticeable 
when  a shallow  source  is  used  and  results  from  the  interference  of  ray  paths 
that  are  upgoing  and  downgoing  at  the  source.)  When  energy  is  allowed 
to  refract  or  reflect  back  into  the  water  column,  the  transmission  loss 
curve  will  show  added  interference  structure.  Hanna's  method  proposes  to 
identify  interference  nulls  in  transmission  loss  CLirves  that  arise  because 
of  refracted  energy  and  to  determine  velocity  and  gradient  information  from 
the  positions  of  these  nulls. 
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VI.  SUPPORT  WORK 


This  chapter  contains  a brief  summary  of  modeling  efforts  conducted 
at  ARL:OT,  undei-  this  contract,  in  support  of  the  bottom  interaction 
acoustic  studies.  This  summary  consists  of  three  sections  devoted  to  (l) 
normal  mode  model,  (2)  an  update  of  the  bottom  loss  model,  and  (3)  an  update 
of  the  ray  trace  model.  All  three  models  discussed  in  this  chapter  were 
developed  at  ARL:Ur  in  the  bottom  interaction  program  to  solve  specific 
bottom  interaction  problems. 

A.  A Normal  Mode  Ifcdel 

A parallel  shooting  technique  together  with  numerical  integration  is 
used  to  solve  the  depth  separated  wave  equation.  The  model  is  restricted 
to  low  frequencies  and  a horizontally  stratified  deep  ocean  but  allows 
arbitrary  sound  speed  profiles  with  a detailed  bottom  description.  Initial 
estimates  for  the  first  eigenvalue  are  obtained  using  harmonic  oscillator 
eigenvalues.  Estimates  for  the  nth  eigenvalues  are  obtained  from  the  ray 
theory  cycle  distance  using  the  (n-l)st  eigenvalue,  and  the  eigenvalue 
spacing  of  the  (n-l)st  and  (n-2)nd  eigenvalues.  The  secant  method  is  then 
used  to  improve  the  eigenvalue  estimates.  Integration  of  the  depth  separat- 
ed equation  is  done  toward  the  axis  from  each  of  the  boundaries  and  is 
matched  at  the  sound  channel  axis,  thus  eliminating  much  of  the  numerical 
Instability  that  is  associated  with  some  modes.  In  addition,  group  veloci- 
ties are  computed  by  evaluating  an  integral  rather  than  by  direct  numerical 
differentiation.  This  method  of  computing  group  velocities  gains  in  both 
speed  and  accuracy  relative  to  differential  methods. 

Auxiliary  programs  compute  and  plot  propagation  loss  (fully  coherent 
and  incoherent  mode  sums)  as  well  as  the  modes  and  sound  speed  profiles, 
ray  and  mode  cycle  distances,  phase  and  group  velocities,  and  other  useful 
quantities.  Software  is  under  development  to  permit  the  use  of  multisensor 
2-dimensional  receiving  arrays. 
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Considerable  attention  has  been  paid  to  producing  a noninteractive 
mode  model  which  is  relatively  fast  and  has  sufficient  accuracy  to  permit 
computation  of  propagation  loss  to  long  ranges.  Execution  time  of  present 
versions  of  this  model  are  approximately  3 to  5 sec  per  mode  on  a CDC  3200 
computer.  This  time  includes  the  time  to  compute  the  mode  at  each  of  the 
several  thousand  depth  points  in  the  mesh.  Each  entire  mode  as  well  as 
eigenvalue,  group  velocity,  mode  attenuation,  etc.,  is  written  on  a magnetic 
tape  for  later  use. 

A version  of  this  model  is  complete  and  is  in  daily  use.  Further 
refinements  are  being  made,  particularly  to  ensure  stability  of  solution 
and  complete  noninteractiveness.  A report  describing  in  detail  the 
characteristics,  advantages,  and  limitations  of  this  model  is  in  prepara- 
tion and  will  be  issued  during  the  next  contract  year  (FY  T?)-^ 

B.  Updates  of  the  Bottom  Loss  Model  (BOTLOSS) 

2 

A previous  report.  Hawker  et  al.,  contains  a brief  description  of  a 

plane  wave  bottom  reflection  loss  model  based  on  nioraerical  solution  of  the 

wave  eq.uation  (program  BOTLOSS ) . Subsequently,  a detailed  technical  report 

2 

on  the  model,  Hawker  and  Foreman,  was  issued  which  described  its  mathemati- 
cal underpinnings,  computational  methods,  and  research  applications.  Since 
the  report  was  released,  BOTLOSS  has  undergone  further  revisions.  Most 
of  these  changes  were  implemented  to  reduce  program  run  time  and  to  provide 
more  flexible  input/output  options.  These  routine  modifications  will  not 
be  described  here  because  they  pertain  only  to  the  details  of  running  the 
program.  However,  there  have  also  been  more  substantive  changes  which 
affect  the  representation  of  the  bottom  sediments  in  the  model.  These 
alterations  are  the  subject  of  this  section. 

1.  Background 

Consider  the  model  ocean  bottom  of  Fig.  VI-1;  the  bottom  consists 
of  horizontally  stratified  Inhomogeneous  fluid  sediment  layers  overlaying 
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an  infinite  homogeneous  substrate  vhich  may  be  a fluid,  a solid,  a perfect 
reflector,  or  a pierfect  absorber.  If  plane,  time  harmonic  acoustic  waves 
inplnge  on  the  bottom  at  a given  grazing  angle  0,  then  the  reflection 
coefficient  is  defined  to  be  the  ratio  of  the  pressure  amplitude  of  the 
upgoing  (reflected)  waves  to  the  amplitude  of  the  duwngoing  waves.  The 
total  press’ire  in  the  water  is  given  by  the  superposition  of  the  upgoing 
and  downgoing  waves.  Assuming  the  water  to  be  homogeneous  and  nonattenuating, 
then  the  time  independent  depth  component  of  pressure,  p^,  may  be  expressed 
as 

iK  z -IK  z 

p(z)  =A(e  ° +Re  where  (l) 

o 0 

z is  depth, 

A^is  the  pressure  amplitude  of  the  downgoing  wave, 

K = cu/c  sin  Q, 

CO  is  2«f, 

is  the  sound  speed  in  water,  and 
R is  the  reflection  coefficient. 


Evaluating  p^and  dp^/dz  at  z=0  gives 

R _ ^ (2) 

^ - 1+S 


6 


P,'(0) 


1 

iK 


(5) 


One  can  therefore  compute  R given  p(o)  and  p'(0).  To  find  these  quantities 
BOTLOSS  must  first  compute  the  pressure  field  numerically  over  the  entire 
interval  from  the  top  of  the  substrate  to  the  water- sediment  interface. 

The  equation  governing  the  pressure  field  in  each  sediment  layer  is 


Pj"(z) 


P.'(z) 


Pj'(z) 


Kj  (z)  Pj(z)  = 0 


(M 


-i6o- 


9 


(5) 


K ^ (z)  =k.^(z)  - cos^  6 
J J ° 


“o'"’  "rr^T*  “j‘"’ 

J 


(6) 


where  a.(z)  Is  attenuation,  p.(z)  is  density,  and  c.(z)  is  sound  speed. 

(]  J J 

At  each  fluid  sediment  layer  interfrce  we  have  the  boundary  conditions 


(7) 


If  the  substrate  is  a fluid  or  solid,  one  can  calculate  the 
pressure  at  the  bottom  of  the  adjacent  fluid  layer  up  to  a constant  using 
given  values  of  substrate  density,  sound  speed,  shear  wave  speed,  and  shear 
and  compressional  attenuations.  The  calculation  is  greatly  simplified  if 
the  substrate  is  a perfect  absorber  or  reflector.  The  undetermined  con- 
stant turns  out  to  be  an  amplitude  coefficient  whose  exact  value  is 
related  only  to  the  pressure  field  normalization.  BOTLOSS  is  not  concerned 
with  the  normalization;  thus  this  coefficient  is  set  arbitrarily  to  unity. 
BOTLDSS  proceeds  to  compute  the  pressure  field  by  solving  the  differential 
Eq.  (4)  with  layer  interface  conditions  (7)  and  (8)  and  initial  values 
given  at  the  sediment -substrate  interface  by  using  a standard  numerical 
routine.  When  the  solution  is  complete,  the  values  of  p^'(o)  and  Pq(0)  are 
known  and  the  program  can  Immediately  calculate  R by  Eqs.  (2)  and  (3). 

2 . Liberalization  of  Sediment  Profile  Descriptions 

Note  that  the  pressure  field  Eq.  (4)  calls  for  the  evaluation  of 

the  profile  functions  p (z),  p '(z),  c (z),  anda  (z).  BOTLOSS  permits  the 

J J J (J 

user  to  select  analytic  profile  functions  appropriate  to  his  needs  for 
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each  layer  from  among  a set  of  available  functions.  The  repertoire  of 
available  functions  has  been  expanded  slightly  and  now  includes 


p.(o)  (constant  density) 


pj(z)  = / (constant  density  gradient) 


p.(o)  e^^  (exponential  density  profile) 
J 


c (z)  = 
J 


c (o)  (constant  sound  speed) 

kJ 


c.(o)  + gz  (constant  sound  speed  gradient) 
J 


c • (o)/ y 1+Pz  (pseudolinear  sound  speed  profile) 

J 

c.(o)/  ^/l+A^(e  '■"'^-l)  (modified  exponential  sound  speed 

profile) 


a (z)  = 
J 


a (o)  (constant  attenuation) 
J 


®.(o)  + gz  (constant  attenuation  gradient) 
J 


a (o)/  (pseudolinear  attenuation  profile) 

J 


a (O)/  v(L+A^(e  ^^-l)  (modified  exponential  attenuation 
“ profile) 

a (O)  c (o)/c  (z)  (soxind  speed  linked  attenuation  profile, 
J J J 


where  A,  p,  and  g are  constant  parameters.  The  modified  exponential  profile 
functions  have  an  extra  parameter  which  causes  the  profile  gradient  to 
decrease  at  great  depths,  in  accordance  with  some  measured  data. 
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These  particular  choices  of  profile  functions  were  made  largely 
because  considerable  research  has  been  done  in  the  past  using  them. 
Typically,  these  functions  were  used  because  they  lead  to  exact  solutions 
or  simple  approximations  in  the  analysis  of  bottom  loss  effects.  It  is 
expected  that  these  functions  will  change  in  the  future,  particularly  in 
response  to  new  Information  about  measiared  profiles. 

It  should  be  noted  that  establishing  the  attenuation  function 
a(z)  in  its  own  right.  Instead  of  linking  the  attenuation  to  the  sound 
speed  profile,  represents  a considerable  departure  from  previous  practice. 
This  development  is  the  subject  of  the  next  section. 

5.  Uncoupling  the  Attenuation  from  the  So\jnd  Speed 

In  most  published  reflection  coefficient  models  and  in  earlier 
versions  of  BOTLOSS,  the  complex  wave  number,  k,  has  been  written  as 


The  attenuation  as  a function  of  depth  is  thus  implicitly  defined  to  be 
inversely  proportional  to  the  sound  speed 

a(z)=a(0)f[^  . (10) 

The  sound  speed  and  attenuation  have  traditionally  been  coupled  in  this 
fashion  because,  for  certain  assmed  forms  of  (real)  c(z),  one  can  find 
analytical  solutions  to  the  wave  equation  and,  by  linking  attenuation  to 
sound  speed  as  in  Eq.  (lO),  one  may  introduce  attenuation  effects  into  the 
solution  without  making  the  analysis  more  difficult. 

Program  BOTLOSS,  however,  solves  the  wave  equation  nvunericallyj 
thus,  there  is  no  need  to  write  k(z)  in  a form  which  renders  the  wave  equation 
analytically  tractable.  The  current  version  of  BOTIDSS,  therefore,  has 
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k(z)  = a3/c(z)  + ia(z) 


(11) 


where  c(z)  and  a(z)  are  independent  of  each  other.  The  nianerical 
integration  algorithm  is  largely  indifferent  to  the  exact  forms  of  c(z) 
and  q;(z)  from  the  standpoint  of  execution  time  and  accuracy}  thus  the  user 
is  free  to  choose  them  to  be  as  physically  realistic  as  possible.  At 
present,  BOTIXDSS  is  equipped  to  handle  constant,  linear,  pseudolinear, 
modified  exponential  attenuation  profiles  and  sound  speed  linked  attenua- 
tion profiles.  This  repertoire  will  almost  certainly  change  as  more  infor- 
mation becomes  available  concerning  actual  variations  of  attenuation  with 
depth . 

Some  preliminary  studies  have  shown  that,  for  certain  types  of  very  thick 
sediments,  the  implicit  attenxiation  function  of  Eq.  (9)  is  inadequate;  this 
results  in  bottom  loss  predictions  which  differ  markedly  from  predictions 
obtained  using  Eq.  (lO)  and  a(z)  chosen  to  fit  measvired  attenuation  profiles 
(see  Chapter  II.  B). 

4.  Reduction  in  Program  Run  Time 


An  aspect  of  the  numerical  Integration  approach  used  in  BOTIOSS 
is  its  relatively  slow  speed,  which  limits  use  of  the  program  to  sediments 
no  thicker  than  about  100  wavelengths  for  reasonable  execution  times.  Most 
of  the  run  time  is  spent  in  the  numerical  solution  of  the  differential 
equation.  Some  of  the  factors  which  influence  integration  time  are  discussed 
in  the  technical  report  on  the  model.  However,  the  difficulty  is  essentially 
that  both  the  pressure  and  its  derivative  are  rapidly  varying  functions  of 
depth  and  therefore  the  numerical  integrator  is  forced  to  take  very  small 
steps  in  order  to  control  errors.  It  was  hoped  that  the  rapid  variation  of 
the  wave  function  could  be  removed  by  writing  the  solution  as  a sum  or 
product  or  other  suitable  combination  of  an  analytical  approximation, 
and,  a slowly  varying  corrector  which  could  be  found  n\jmerically.  For 
example,  one  might  write 
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p(z)  = PjwKB^^^ 

where  p(z)  is  the  desired  exact  solution,  Pjyj^C^)  is  the  JWKB  approximation 
to  the  solution,  and  T(z)  is  a numerically  determined  corrector.  Introduc- 
ing this  expression  into  the  wave  equation  produces  a new,  nonlinear 
differential  equation  with  solution  T,  which  is  (it  is  hoped)  slowly 
varying.  Several  such  schemes  were  examined  but,  so  far,  these  attempts 
have  been  only  partially  successful.  Typically,  whenever  the  behavior  of 
p(z)  was  approximately  exponential,  the  numerical  integration  was  greatly 
accelerated  because  the  Integrator  could  take  much  larger  steps;  however, 
the  problems  of  integrating  in  oscillatory  portions  of  the  field  remained. 
Roughly  speaking,  the  reason  for  the  failures  in  the  oscillatory  regions 
was  that  the  approximate  solution  contained  errors  in  both  the  magnitude 
and  phase  of  oscillation;  although  these  errors  were  usually  small,  the 
corrector  could  not  compensate  for  both  of  them  simultaneously.  These 
schemes  based  on  the  JWKB  approximation  were  abandoned  after  a short  in- 
vestigation but  they  may  merit  a closer  look,  especially  since  any  success- 
ful techniques  which  are  developed  will  likely  be  applicable  to  a much 
broader  range  of  problems. 

A modification  of  this  approach  which  holds  considerable  promise 

is  to  use  in  place  of  the  JWKB  term  an  asymoptotic  approximation  due  to 

4 

Langer.  The  utility  of  this  approach  in  acoustic  normal  mode  calculations 

5 

has  already  been  investigated  at  this  laboratory  by  Mitchell.  It  is  anti- 
cipated that,  when  a clear  need  arises  for  significantly  decreasing  the 
execution  time  of  the  program,  this  method  will  be  implemented. 

C.  Developments  in  Acoustic  Ray  Theory  Programs 

2 

The  final  report  for  1975  on  the  bottom  interaction  study,  sponsored 
by  NAVEIEX,  Code  520,  includes  a description  of  a computer  program  devel- 
oped at  ARL:UT;  this  program  numerically  locates  eigenrays  (rays  connecting 
a sotirce  and  receiver)  and  prints  their  travel  times,  laianch  angles,  bottom 
reflection  angles,  and  other  ray  path  Information.  This  program,  RANGER  BP, 
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has  a limited  capacity  to  find  eigenrays  which  travel  through  ocean  bottom 
sediments.  By  providing  a detailed  path  history  of  each  eigenray  arriving 
at  a given  receiver,  RANGER  BP  has  aided  in  understanding  and  predicting 
the  effects  of  multipath  arrivals.  Its  capacity  to  trace  rays  through 
the  bottom  made  it  useful  in  the  investigation  of  acoustic  techniques  for 
probing  the  bottom  to  determine  such  sediment  properties  as  sound  speed  and 
layer  thickness. 

Since  1975^  ray  theory  modeling  capability  at  ARL;OT  has  been  upgraded 
and  expanded  to  include  several  programs,  collectively  referred  to  as 
RANGER.  The  most  important  development  has  been  the  added  ability  to  com- 
pute the  acoustic  pressvire  contribution  of  each  eigenray  to  the  total 
pressure  at  a receiver. 

All  programs  in  the  RANGER  package  assiame  a range  invariant,  horizontally 
stratified  ocean  acoustic  environment.  The  sound  speed  profile  is  approxi- 
mated by  piecewise  linear  segments,  for  which  analytical  expressions  are 
available  for  computing  travel  time,  path  length,  horizontal  distance  (range), 
and  range  derivatives  as  functions  of  the  ray  "launch  angle."  These  geo- 
metric ray  path  quantities  permit  calculation  of  the  ray  pressure  at  the 
receiver  within  the  limits  of  the  ray  theory  approximations. 

The  computational  methods  used  in  RANGER  are  discussed  in  detail  in  a 
separate  report  (Foreman  ) . The  descriptions  of  the  programs  currently  in 
the  RANGER  package  are  summarized  below  with  emphasis  on  their  applications 
to  the  bottom  interaction  study. 

RANGER  BL  - Computes  and  plots  propagation  loss  as  a function  of  range. 
Intensities  are  calculated  by  locating  eigenrays  and  summing  their  indivi- 
dual contributions  until  the  total  intensity  begins  to  converge,  whereupon 
the  search  for  additional  eigenrays  is  discontinued.  Both  incoherent  and 
fully  coherent  intensity  summations  are  performed.  Tables  of  bottom  loss 
and  phase  as  functions  of  grazing  angle  are  used  to  account  for  bottom 
interactions.  (No  rays  are  traced  through  the  sediment  in  this  version.) 
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The  ability  to  fully  account  for  phase  interference  effects  has  been  crucial 
in  accounting  for  some  effects  observed  with  receivers  in  close  proximity 
to  each  other  or  to  the  ocean  surface  or  bottom. 


RANGER  MP  - Prints  an  account  of  the  ray  path  history  of  each  eigenray 
which  arrives  at  the  receiver  within  a specified  time  interval  of  the 
earliest  arrival.  It  predicts  multipath  arrival  times,  pressure  levels, 
arrival  angles,  tin-ning  depths  or  surface/bottom  reflection  angles,  phase 
shifts  due  to  ray  interactions  with  the  siorface  or  bottom  or  caustics, 
and  a general  description  of  the  ray  (upgoing/downgoing  at  the  source/re- 
ceiver, number  of  shallow  and  deep  turning  points,  etc.).  This  program 
is  particularly  useful  for  Identifying  the  ray  paths  which  determine  the 
multipath  arrival  structure  in  experimental  data. 

RANGER  BP  - Prints  ray  histories  similar  to  RANGER  MP.  It  has  a 
limited  capacity  to  trace  rays  through  sediment  layers.  It  can  be  used 
to  predict  multipath  arrival  times  for  rays  which  undergo  only  one  bottom 
interaction.  RANGER  BP  does  not  perform  intensity  calculations. 

JiAYEAM  - Plots  ray  paths  through  the  ocean  for  given  launch  angles. 

RTHETA  - Plots  ray  range  as  a function  of  launch  angle  for  a given 
number  of  deep  turning  points  (bottom  reflections  or  deep  refractions). 
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Code  2627 


PREDEDINO  PAGS  BULtBUNOT  FlUiSi) 


L.. 


173 


Dist.  List  for  ARL-TR-TT-2T  under  Cont.  N00039-76-C-0081,  Cont'd 


Copy  No. 

22-21+ 

25 

26 

27 

Commander 

Naval  Ocean  Systems  Center 
Department  of  the  Navy 

San  Diego,  CA  92152 

Attn:  E.  Tunstall  (Code  531) 

R.  R.  Gardner  (Code  53) 

M.  A.  Pederson  (Code  307) 

H.  P.  Bucker  (Code  5311) 

28 

29 

Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 

Washington,  DC  20362 

Attn:  C.  D.  Smith  (Code  06HI) 

A.  P.  Franceschetti 

30 

Chief  of  Naval  Operations 
Department  of  the  Navy 

Washington,  DC  20350 

Attn:  R.  S.  Winokur  (OP95E1) 

31 

Chief  of  Naval  Material 

Department  of  the  Navy 

Washington,  DC  20360 

Attn:  G.  R.  Spalding  (Code  08T2lt) 

32 

Commander 

Naval  Intelligence  Support  Center 
1+301  S\iitland  Road 

Washington,  DC  20390 

Attn : Code  222 

33 

Oceanographer  of  the  Navy 

Hoffman  Building 

200  Stovall  Street 

Alexandria,  VA  22332 

Attn:  CDR  E.  Young 

3U  Commander 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

35  Commander 

David  M.  Taylor  Naval  Ship  Research  and  Development  Center 
Bethesda,  MD  2003^ 
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36 

37 

38 

39 


ho 

Ul 

1*2 


1*3 

1*1* 

1*5 

U6 


Naval  Oceanographic  Office 
Department  of  the  Navy 
Washington,  DC  20373 
Attn:  W.  H.  Geddes 

Commander 

Naval  Air  Development  Center 
Department  of  the  Navy 
Warminster,  PA  l897l* 

Attn:  P.  Van  Schuyler  (Code  2052) 
C.  L.  Bartherger 
P . Haas 


Commander 

New  London  Laboratory 
Naval  Underwater  Systems  Center 
Department  of  the  Navy 
New  London,  CT  06320 
Attn:  F.  R.  DiNapoli 

R.  Martin 

Assistant  Director 
Ocean  Control  DDR&E 
Room  3D,  10 1*8  Pentagon 
Washington,  DC  20301 

OASN  (R,E  & S) 

Room  1*D,  71*5  Pentagon 
Washington,  DC  20301 
Attn:  G.  A.  Cann 

Advance  Research  Project  Agency 
Acoustic  Research  Center 
Moffett  Field,  CA  91*035 
Attn:  T.  Kooij 

Superintendent 
Naval  Postgraduate  School 
Monterey,  CA  939l*0 
Attn:  Library 

Commanding  Officer 

Naval  Air  Systems  Command 

Department  of  the  Navy 

Washington,  DC  20360 

Attn:  CDR  J.  Messegee  (Code  PMA-261*) 


175 


Dist.  List  for  ARL-TR-TT-2T  'jnder  Cont . N000  39-T6-C-0081  , Cont'd 


Copy  No. 
UT-U8 


1*9 


50 

51 


52 

53 


51* 


55 


56 


57 

58 


Commanding  Officer  and  Director 
Defense  Documentation  Center 
Defense  Services  Administration 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria,  VA  22314 

Arthur  D.  Little,  Inc. 

Acron  Park 
Cambridge,  MA  02140 
Attn:  G.  Raisbeck 

Woods  Hole  Oceanographic  Institution 
86-96  Water  Street 
Woods  Hole,  MA  02543 
Attn:  E.  E.  Hays 
G.  Frisk 

Science  Applications , Inc. 

8400  Westpark  Drive 
McLean,  VA  22101 
Attn : J . Hanna 
C.  Spofford 

Pennsylvania  State  University 
Applied  Research  Laboratoiy 
P.  0.  Box  30 

State  College,  PA  168OI 
Attn:  S.  McDaniel 

Underwater  Systems,  Inc. 

World  Building 
8121  Georgia  Avenue 
Silver  Spring,  MD  20910 
Attn:  M.  S.  Weinstein 

University  of  California,  San  Diego 

Marine  Physical  Laboratory  of  the  Scripps  Institution  of 
Oceanography 
San  Diego,  CA  92132 
Attn:  V.  Anderson 

Tracer,  Inc. 

1601  Research  Boulevard 
Rockville,  MD  20850 
Attn:  J.  Gottwald 

A.  Wittenbom 
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59 

60 


61 

62 

63 

6h 

65 


66 

67 


Bell  Telephone  Laboratories,  Inc. 
Whippany  Road 
Whippany,  NJ  07981 
Attn:  S.  A.  Kramer 

Planning  Systeirn  Incorporated 
7900  Westpark  Drive,  Suite  507 
McLean,  VA  22101 
Attn:  L.  Solomon 

TRW,  Inc. 

TRW  Defense  & Space  Systems  Group 

Washington  Operations 

7600  Colshire  Drive 

McLean,  VA  22101 

Attn:  R.  T.  Brown 

Xonics,  Inc. 

6837  Hayvenhurst  Avenue 
Van  Nuys,  CA  91^*06 
Attn:  J.  Wilson 

SUT'RON  Corporation 
Suite  700 

1923  North  Lynn  Street 
Arlington,  VA  22209 
Attn : C . Dabney 

Daubin  Systems  Corporation 
lOL  Crandon  Boulevard 
Key  Biscayne,  FL  331*^9 
Attn : S . Daubin 

A.  0.  Williams,  Jr. 

Dept,  of  Physics 
Brown  University 
Providence,  RI  02912 

Defense  Scientific  Establishment 
HMNZ  Dockyard 
Devonport,  Auckland 
New  Zealand 
Attn:  K.  M.  Guthrie 
R.  N.  Denham 
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68 

69 

70 

71 

72 

73 

7lt 

75 

76 

77 

78 

79 

80 

81 
82 
83 
QU 


University  of  Auckland 
Physics  Department 
Auckland,  New  Zealand 
Attn:  A.  C.  Kibhlewhite 
G.  Bold 

The  Catholic  University  of  America 
6220  Michigan  Avenue  NE 
Washington,  DC  20017 
Attn:  H.  M.  Uberall 

TRACOR,  Inc. 

1601  Research  Boulevard 
Rockville,  MD  20850 
Attn:  R.  J.  Urick 

Geophysical  and  Polar  Research  Center 
Department  of  Geology  and  Geophysics 
The  University  of  Wisconsin 
Madison,  WI  53700 
Attn:  C.  S.  Clay 

D.  C.  Stickler 
Courant  Institute 
251  Mercer  Street 
New  York,  NY  10012 

Glen  E.  Ellis,  ARL:UT 

Terry  L.  Foreman,  ARL:UT 

Ruth  Gonzalez,  ARL:UT 

loyd  D.  Hampton,  ARL:UT 

Kenneth  E.  Hawker,  ARL:UT 

John  Lemmon,  ARL:UT 

Stephen  K.  Mitchell,  ARL:UT 

Susan  Payne,  ARL:UT 

Steve  Rutherford,  ARL:UT 

Jack  A.  Shooter,  ARL:UT 

Chris  Tindle,  ARL:UT 
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85 

Donna  Venable,  ARL:UT 

86 

Paul  Vidmar,  ARL:UT 

87 

Reuben  H.  Wallace,  ARL:UT 

88 

Steven  L.  Watkins,  ARL:UT 

89 

Library,  ARL:UT 

90-110 

Reserve,  ARL;UT 
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